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§ ' Abstract 

^ : 

^ ■ We have presented a detailed study of twist-2 and twist-3 hght-cone distribution amphtudes 
^ ■ of l^Pi and l^Pi axial-vector mesons, based on QCD conformal partial wave expansion. 

Applying equations of motion, the twist-three two-parton light-cone distribution amplitudes 
p i' can be expressed in terms of leading-twist and twist-three three-parton light-cone distribu- 
p tion amplitudes. The relevant G-parity invariant and violating parameters, containing the 



corrections due to the SU(3) breaking effects, are evaluated from the QCD sum rule method. 
The results for axial- vector decay constants of l^Pi states are presented. The values of tensor 
CN ■ decay constants and Gegenbauer moments of the leading twist distribution amplitudes for 
^ ■ l^Pi states are updated. Using Gell-Mann-Okubo mass formula, the mixing angle for the fg 

On . and /i of l^Pi states is 03 ~ 38°, and that for and hi of l^Pi states is ftp^ ~ 10°. The 
■ detailed properties for physical states /i (1285), /i(1420), 170), and /;,i(1380) are given, 
l/-^ ■ Assuming the mixing angle between Kia and Kib to be 9k = 45° or —45°, we also give the 
detailed study for ii"i(1270) and i^i(1400). Using the conformal partial wave expansion, we 
O ■ obtain the models for light-cone distribution amplitudes, containing contributions up to con- 
formal spin 9/2. It is interesting to note that some distribution amplitudes have significant 
^ . asymmetric behaviors, which should be phenomenologically attractive. 
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1 Introduction 



For an energetic light hadron moving nearly on the light-cone in the reference frame, its dis- 
tribution amplitudes can described by a set of light-cone distribution amplitudes (LCDAs) 
[U E]. The LCDAs are governed by the special collinear subgroup SL{2,'R) of the confor- 
mal grouplH [2]. There are four generators for the 5'L(2,M) group. Three of generators 
describe the generalized rotations in the SL{2,M.) space, where the rotational invariance is 
characterized by the so-called conformal spin j, in analogy to the orbital quantum number 
in quantum mechanics of having spherically symmetric potential. The remaining generator 
counting the collinear twist commutes with the rest generators in the SL(2, M) group; thus, 
in the mathematical spirit, the concept of "collinear twist" is equivalent to the "energy" in 
quantum mechanics, and, in other words, the SL{2, M) group are satisfied by the same Lie 
algebra as the 0(3) group. 

The conformal partial wave expansion of a light-cone distribution amplitude is fully analo- 
gous to the partial wave expansion of a wave function in quantum mechanics. Each conformal 
partial wave is labeled by the specific conformal spin j. On the other hand, for a hadron mov- 
ing nearly on the light-cone, the transverse separation of the partons corresponding to the 
two-dimensional transverse plane (0, Xi, X2, 0) is characterized by the SL{2, C) group. More- 
over, the collinear SL{2, M) and transverse SL{2, C) groups are not independent because they 
share the same dilation generator. Integrating out the transverse degrees of freedom, which 
is governed by renormalization group (RG) equation, yields the scale-dependent behavior 
for each conformal partial wave amplitude of a hadronic LCDA. The conformal invariance 
in QCD exhibits that for leading twist LCDAs there is no mixing among conformal partial 
wave amplitudes with different conformal spins to leading logarithmic accuracy, while, for 
higher twist LCDAs, because there may exist multiple operators with the same conformal 
spin, the corresponding conformal partial amplitudes with the same conformal spin can mix 
as the scale is changed. Roughly speaking, to leading logarithmic accuracy, the anomalous 
dimensions rise logarithmically with conformal spin. This implies that contributions of con- 
formal operators with high conformal spins are suppressed at large scale; for instance, the 
well-known leading twist LCDA of the pion reads — oo) 6u{l — u), i.e., the 

asymptotic amplitude with the lowest conformal spin 2, where u is the momentum fraction 
carried by the quark (or anti-quark) in the pion. It is interesting to note that in the BFKL 
approach [31 IH El E], the QCD scattering amplitude in the large energy limit is related to 
the 5'L(2,C) group (see, e.g., the discussion in Ref. [2]). 

Understanding the hadronic LCDAs is very important since in the QCD description 
of various exclusive processes the amplitudes can be represented as the hadronic LCDAs 
convoluted with the interacting hard kernel. The studies for the hadronic LCDAs were first 
done by Chernyak and Zhinitsky [7j. They used the QCD sum rule approach to calculate 
the moments of LCDAs of vr, K, p, K* and some baryons. Ten years later, a conformal 
description for multi-parton LCDAs was presented in more details in Ref. [Tj. Most recently, 
using conformal expansion together with equations of motion (EOMs), the different twist 
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LCDAs for the light pseudoscalar and vector mesons have been systematically studied in 

Refs. iiiniio]. 

In the present paper, we devote to examining twist-2 and twist-3 LCDAs of light axial- 
vector mesons with quantum numbers l^Pi, l^Pi and their mixtures. The motivation of 
this study is as follows. LCDAs of axial-vector mesons should be important for (exclu- 
sive) phenomenologies involving axial- vector mesons. 5° 7^1(1270)7 has recently been 
measured by Belle [11], whereas the real physical states -ft'i(1270) and -ft'i(1400) are the 
mixtures of ideal l^Pi (Kia) and l^Pi (Kib) states. The first charmless hadronic B decay 
involving a l^Pi meson that has been observed is P° a^(1260)7r^ [12l[l3], which is rele- 
vant to measurements of the unitarity triangle a = 02 of the Cabibbo-Kobayashi-Maskawa 
(CKM). BaBar has recently extracted the effective value aes from the measurements of CP- 
violating asymmetries in the decay P° — ^ a]''(1260)7r''' [13], where the bound on the difference 
Aa = a — a^s can be constrained by using the broken SU(3) flavor symmetry [15]. On the 
other hand, B{B ^ V A, A A) can be used to test QCD annihilation topology in B decays 
and probe the new-physics signals, where V = vector meson and A = axial-vector meson 
[MlIITj. 

So far, there is no literature for the calculations of LCDAs of axial-vector mesons except 
my previous work about leading-twist LCDAs of l^Pi states [18j. In Ref. [18], we have 
calculated the Gegenbauer moments (with conformal spins 0, 1, and 2) of the leading twist 
LCDAs for the 6i(1235), /ii(1380) and Kib as well as their tensor decay constants. In this 
paper, we will derive all parameters, relevant to the two- and three-parton LCDAs of twist- 
2 and -3, with conformal spin up to 9/2. We take into account SU(3) breaking effects in 
the sum rule calculations for G-parity invariant and violating parameters, where the latter 
parameters should vanish in the SU(3) limit. The decay constants for axial- vector mesons 
can be obtained by using the sum rules computed in Ref. [19], where the authors focused 
on the studies of the masses and used a small value of as- Considering the current value of 
as and the renormalization group improvement for sum rules, we also examine the masses 
of the complete l^Pi set, which may offer information for the quality of the sum rules. Our 
results for parameters relevant to LCDAs of axial-vector mesons are essentially original. To 
determine the relative signs for parameters, one of the interpolating currents in the two- 
point correlation function is adopted to be the local axial-vector current (or the local tensor 
current) for I'^Pi states (or l^Pi states). However, because the resultant f^sp^ sum rule does 
not show reliable quality, we thus resort to the diagonal correlation function to calculate the 
sum rule for this parameter. Fortunately, the related calculations can be found in Ref. |2Uj . 
where the sum rule is instead used to analyze the couplings for pseudoscalar mesons. The 
detailed results are given in Subsec. I5.1.3[ Except the formula f^ip^ can be adopted directly 
from Ref. [20], we do not find any explicit results can be applied directly, but only two 
equations can be compared with the results in the massless quark limit in the literature. See 
discussions before and after Eq. (11211) and after Eq. (12441) . 

We also update the numerical results for the decay constants and leading twist LCDAs of 
the l^Pi states, due to the following reasons. First, here we consider the uncertainties of the 
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condensates and strange quark mass, which were not fully included in Ref . [TS] . Meanwhile, 
we re-examine the stability of the sum rules with requirement that the contributions from 
excited states and from the highest dimension term in OPE expansion can be well under 
control within the working Borel window. Second, in Ref. [18] we assumed that the prop- 
erty of /ii(1170) is the same as that of 6i(1235), and the quark content of /;,i(1370) is ss. 
Instead, we study the pure l^Pi states here. We estimate the singlet-octet mixing angle by 
means of Gell-Mann-Okubo mass formula (see the detailed discussions in subsection 14. ip . 
Consequently, the physical properties for the real states /ii(1170) and /ii(1380) are obtained. 
Third, the G-parity violating Gegenbauer moments were not evaluated correctly in Ref. 
because the corrections arising from the l^Pi states had been ignored. In the sum rule 
calculations, for G-parity violating parameters, the contributions originating from l^Pi and 
l^Pi are always of the same order. 

In the quark model, a l^Pi (1^-Pi) meson is represented as a constituent quark- ant iquark 
pair with total spin S* = 1 (0) and angular momentum L = 1 (1). Nevertheless, a real 
hadron in QCD language should be described in terms of a set of Fock states for which each 
state has the same quantum number as the hadron, and the leading-twist LCDAs are thus 
interpreted as amplitudes of finding the meson in states with a minimum number of partons. 
Interestingly, due to the G-parity, the leading-twist LCDA $^ ($||) of a l^Pi (1^-Pi) meson 
defined by the nonlocal tensor current (nonlocal axial-vector current) is antisymmetric under 
the exchange of quark and anti-quark momentum fractions in the SU(3) limit, whereas the 
$11 ($±) defined by the nonlocal axial- vector current (nonlocal tensor current) are symmetric. 
The large magnitude of the first Gegenbauer moment of $|| ($_l) could have a large impact 
on the longitudinal fraction of factorization-suppressed B decays involving a l^Pi (or l^Pi) 
meson evaluated in the QCD factorization framework [161 E]- Furthermore, $^ is relevant 
not only for exploring the tensor-type new-physics effects in B decays [T6| but also for 
B Ki-f studies 

This paper is organized as follows. A brief introduction for conformal partial wave ex- 
pansions of LCDAs is provided in Sec. [H where we first introduce the SL(2, M) group, and 
then discuss the asymptotic LCDAs, two-parton and three-parton distribution amplitudes. 
In earlier days, Appell's polynomials were used to form the conformal basis in expanding 
three-parton distribution amplitudes. However, the Appell's polynomials are not mutual 
orthogonal (see, e.g., Refs. [21 [22]). Following Ref. [21 [221 123], we thus use the orthogonal 
basis, which will be introduced in Sec. 12.21 to expand the three-parton LCDAs of twist-3. 
The detailed coUinear twist properties of LCDAs are collected in Appendix [XI 

Sec. [3] presents the notations, that we use in the paper, and definitions of LCDAs of 
axial-vector mesons. Moreover, we apply the EOMs to obtain the relations among the 
twist-two and twist-three light-cone distribution amplitudes, so that we can use a minimal 
number of independent nonperturbative parameters to describe the distribution amplitudes. 
EOMs in the QCD perturbative theory respect all symmetries given in the classical level; 
in other words, conformal operators related by EOMs should have the same conformal spin. 
Appendix [B] contains some relevant formulas of EOMs for the present work. 
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Using the QCD sum rule technique [23], Sees. [Hand are devoted to the determination of 
relevant parameters of leading-twist and twist-3 three-parton LCDAs, respectively. In Sec. HI 
the physical properties of the axial-vector mesons such as the quark contents and decay 
constants are discussed. The relevant inputs in the calculation are collected in Appendix O 
Sec. [6] contains explicit models for the twist-2 and twist-3 LCDAs of the axial- vector mesons. 
Finally, we summarize in Sec. [71 

2 Conformal Partial Wave Expansions of LCDAs 

The properties of fields living on the light-cone satisfy the SL{2, M) group, which is a coUinear 
subgroup of the full conformal group. The SL{2,M.) group is governed by four generators 
P+, M_+, D and K_. Here and below we introduce two light-like vectors and which 
satisfy n^nf^ = n^n^ = and n^n'^ = 1, so that for a general vector A we define A- = A^n^ 
and = A^n^. In analogy to the space rotation in quantum mechanics, the three linear 
combinations of the S'L(2,M) generators can describe hyperbolic rotations, and have the 
following commutation relations: 

[Lo,L^] = TL^, [L_,L+] = -2Lo, (1) 

where 

L+ = -?PX, L_ = (V2)KX, Lo = -(V2)(D-M^,n'^n^), (2) 

with P^, K^, D and "M.^^ being the translation, special conformal transformation, dilation, 
and Lorentz generators, respectively [25j. Following the standard quantum mechanical tech- 
nique, the above generators acting on a field which lives along the light-cone (in 
other words, the state depends only on = x_n^), yield 

[L+, <|.(x_)] = ^ , (3) 

[L_, <l>(a;_)] = {x'id+ + 2jx_) ^{x^) = L_^x^) , (4) 

[Lo, = - + j) <l>(a;_) = Lo$(a;_) , (5) 

where j = [i + s)/2 with i and s being the field's scaling dimension and spin projection on 
the light-cone, respectively, j is called the conformal spin since 

[L2, = j(j - imx^) = (6) 

where 

= + + = - Lo + L+L_ . (7) 
The remaining generator of the SL{2, M) is symbolized as E = i(D + M__|_) and satisfies 

[E,^x^)] = (8) 

[E,L2]=0, [E,Lo]=0, (9) 
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where i — s = t is called the collinear twis^ of the field <l'(x_). See Appendix Rl for further 
discussions. Therefore the role of E is analogous to the Hamiltonian in quantum mechanics, 
and the twist corresponds to the eigenenergy of the Hamiltonian. In short, a LCDA defined 
by (non-local composite) field operators of a given collinear twist can be expanded by using 
the basis, which is made of the eigenstates of and Lq, as one does the partial wave 
expansion in quantum mechanics. Such expansions for LCDAs are so-called conformal partial 
wave expansions. In QCD, in the large scale limit and to leading logarithmic accuracy, we 
have i = i^'^^ + 7, where i^^^ and 7 are the canonical and anomalous dimensions of the field, 
respectively. 

A generic multi-parton distribution amplitude is defined through the matrix element of 
a multi-local operator 

mU^rn^) ■ ■ ■ ^i{xi^)\p) ~ f [du] e-f(-i«i+-+-'"-™)</,(ui, ■ ■ ■ ,M^) , (10) 



where J[du] = J^dui ■ ■ -J^dum Uj — 1), = are momentum fractions carried by 

partons "i", and the quantum fields $j(xj_) five on the light-cone. (p{ui) can be expanded in 
the Hilbert space by using the basis defined by and Lq. In the expansion, the term with 
the lowest conformal spin and non-zero coefficient is the so-called the asymptotic distribution 
amplitude, which is usually dominant at the large scale for the general case [HIS], 

r(2ji) • ■■T{2jm) 
where the normalization is chosen such that J[dui] 0as(wi) = 1. 

2.1 Two-parton distribution amplitudes 

The study for conformal expansion of two-parton distribution amplitudes has been given by 
Ohrndorf p6] (see, e.g., the discussion in Ref. |2|). Consider 

(0|<l>f (X2)$f (xi)b) = [ [(iM]e-*^'("i"i+"^"^V(wi,W2;/i). (12) 



Using the conformal basis, the two-parton distribution amplitude with a given collinear twist 
can be written as 

00 

0(Mi,M2;/^) = 0as 0j-ii-i2(/^)Pf-ii-i2(^l'"2) • (13) 



j-jl-j2=0 



where 



N 



"1 V "''ij \ XX 

P--(«„«2)= E ^,-,r( li-^ = ( ,(14) 



^For simplicity, it will be just denoted as "twist" in this paper. Conventionally, we call £'^^'^ — s as the 
twist, instead of using i — s, where f'^™ is the canonical dimension. 
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with P'^^[x) being the so-called Jacobi polynomials. Here the conformal spins of ^■j^ and 
$2^ are ji and j^, respectively. 

As for the case of the twist-2 (leading-twist) distribution amplitude of an axial-vector 
meson, we have the conformal spins ji = j2 = 1 for both quark and anti-quark fields. The 
leading-twist LCDA thus reads 

oo 

(j){ui,U2;n) = QuiU2^ai{n)Cf^'^{ui - U2), (15) 
1=0 

where Pj^2i'^i ~ ""2) ~ Cf^'^{ui — U2) with j — 2 = I. 



2.2 Three-parton distribution amplitudes 

We follow the method suggested by Braun et al. |22j, who first applied it to baryon distri- 
bution amphtudes, to construct the three-parton distribution amplitudes of an axial-vector 
meson. The traditional choice is to expand three-parton distribution amplitudes in terms 
of Appell's polynomials. However, this basis is inconvenient for calculations since Appell's 
polynomials are not mutually orthogonal. An orthonormal "conformal basis" of functions 
J^^n^{ui) = {—l)^Yj\j'^^^{ui)/2 can be defined by requiring that the total three-parton con- 
formal spin J = ji + j2 + j3 + N of an eigenstate is fixed, and, moreover, the definite 
value of the conformal spin of the two-parton channel (12) is given as j = ji + j2 + n with 
n = 0, ■ ■ ■ , N . The results of Yj^'^^^ can be written as 

Yjf\u,) = (1 - n3)---- PS^Z^f'-'^il - 2%) Pj^^I^f (f^) • 

Note that this basis of Y'jj^''^(uj) can be related to the basis of Yjj'^^\ui) (or Yjj^^\ui)) 
through the Racah 6j-symbols of SL{2, M) group [221 12] • 

For axial-vector mesons, the three-parton LCDAs correspond to the higher Fock states 
consist of the quark, anti-quark, and gluon. We refer to ji = 1, j2 = 1 and js = 3/2 as the 
conformal spins of the quark, anti-quark and gluon, respectively, in an axial-vector meson. 
Using Eq. ([TB]) . we redefine the orthonormal basis functions, 

^'l^'iu,) = ^-^Y^'^'(u^) = tlUn _ Pl?'l+^")(1 - 2u;) Pi^-^) ( ^^i^ 



J'Nn"i^^d = ^Yy/'-\u,) = ^(1 - P^N-T^l - 2U,) P^^ [f^) . (17) 

where J is the total three-parton conformal spin, j is the total conformal spin of the quark 
and anti-quark, N = J — 7/2 and n = j — 2 with N,n = 0, 1, ■ ■ ■, and N > n. Thus, the 
meson's three-parton distribution amplitudes of twist-three can be represented as 

00 N 
Ar=0 n=0 
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If one takes into account the conformal spins of the distribution amphtudes up to order of 
J = 9/2, the twist-3 LCDA readfl 



(22) 



3 Definitions of light-cone distribution amplitudes 

In what follows, we define z = y — x with z"^ = 0, and introduce the light-like vector 
Pfi = Pii~TT^\z^j_/ {2P-z) with the axial- vector meson's momentum = m\. The polarization 
vector e^^^ of the axial-vector meson can be decomposed into longitudinal and transverse 
projections as [9l [10] 



II M 



(A) 



p ■ Z 



2p ■ z 



(23) 



respectively. In QCD description of hard processes involving axial vector mesons, one en- 
counter bilocal operators sandwiched between the vacuum and the meson. 



{AiP,\)My)T[y,xMxm, 



(24) 



where F is a generic notation for the Dirac matrix structure and the path-ordered gauge 
factor is 



[y,x] = Pexp 



igj dt(x-y)^A''{tx + {l-t)y) 



(25) 



This factor is equal to unity in the light-cone gauge which is equivalent to the fixed- 
point gauge (x — y)^A^{x — y) =0, which is also called the Fock-Schwinger gauge, as 
the quark-antiquark pair is at the light-like separation. For simplicity, here and below we 
do not show the gauge factor. For fi and /ii, the operators in Eq. fl2^ correspond to 



^The contribution of conformal spin 11/2 to the meson's three-parton distribution amplitudes of twist- 
three is 



360uiU2Ui 



(I8M3 - I6U3 + 3)^2,0 + 3(1 - 3m3)(ui - U2)uji.i + ^[(l ^ ""3)^ - 5uiU2]wo,2 



Comparing the above result with that given by Braun and Filyanov [T] , we thus find 



2 3p 1 BF 



° BF ^ BF 



(19) 

(20) 
(21) 



Note that for uji^i there is no corresponding term given in Ref. [T]. 
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{u{y)Tu{x) + d{y)Td{x) + s{y)Ts{x))/ ^/3, while for /§ and hg, the relevant forms of oper- 
ators are {u{y)Tu{x) + d(y)Td{x) — 2s{y)Ts{x)) /^/6. In the present study, we adopt the 
conventions = da + tgsA''^X^/2, G^p = (l/2)e„;3^,G'^^ e^^^s = _i. 

In general, the LCDAs are scheme- and scale-dependent. One can catalog the distribution 
amplitudes into two classes: (i) chiral-even LCDAs for which the relevant current operators 
involve chirality- conserving structures F = {7^,7^75}, and (ii) chiral-odd LCDAs for which 
the operators contain chirality- violating structures T = {o"^^(75), 1}. 



3.1 Two-parton distribution amplitudes 

The chiral-even LCDAs are given by[f| 

{A{P,\)My)j,j,q,{x)\0) =ifAmA / due^^^^y^'^^^^Hp, $|| (n) + e^^^^ (w) 

1 e*^^)^ "1 

^^m7:3^"^a^3(m) \ , (26) 



2 ^{pzy 

{A{P,\)\q,{y)lM^m =-^fAmAt,.p.tll)fz'' /" e^("^^+"^^') ^i^, (27) 



where u and u = 1 — u are the momentum fractions of qi and q2 in the axial- vector meson, 
respectively. The chiral-odd LCDAs are defined by@ 



+ ^^^y KVii^v - VvZ^) (m) 



, 1 / *(A) *(A) X ^\ L / \1 



(2^ 



{A{P, X)My)l,q2{xm = /im^(6*Wz) d«e^("™+^?'-) . (29) 

Here $||, $^ are of twist-2, g^_^\ g^^\ h^^\ h^^^ of twist-3, and (73, of twist-4. In SU(3) limit, 
due to G-parity, 5'^'', and are symmetric [antisymmetric] under the replacement 

M — > 1 — M for the V^Pi [1^-Pi] states, whereas h^\^\ and /is are antisymmetric 

[symmetric]. In other words, in the SU(3) limit it follows that 

1 pi i>i i>i 

du^±{u)= / duh\l\u)= / duh\l'\u)= / duh-i{u) = (30) 
Jo Jo Jo 



g]^' and g\^' given in Ref. [TB] is respectively redefined to be and g]^' in this paper, 
/ill*' given in Ref. [18] is redefined to be h^^^^ in this paper. 
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for l^Pi states, but becomes 

[ = [ dugf{u)= [ dugf{u)= [ dugaiu) = (31) 

Jo Jo Jo Jo 

for l^Pi states. The above integrals are not zero if m^j 7^ rug^, and the detailed results are 
summarized in Sees. 16.11 and 16. 2[ We will not further discuss the twist-4 LCDAs, and 
/13, below. For convenience, we therefore normalize the distribution amplitudes of the l^Pi 
[l^Pi] states to be subject to 



(32) 



/ du^\\{u) = 1 / du^j_{u) = 1 
Jo I Jo 

and take /ap^ = /ap^ [/ip^ = f^^{ii = 1 GeV)] in the study, such that we have 

(l3Pi(P, \)MOK.lMOm = fs\4'''''' (ef - e:Wp,), (33) 
(liPi(P, A)|gi(0)7^75g2(0)|0) = ^/ip^aj''^^ mip, ) , (34) 

where Oq ' and Oq are the Gegenbauer zeroth moments, defined in Eqs. fl60|) and fj48|) . 
and vanish in the SU(3) limit. Using the results given in Eqs. ( fTTI) . ( fT2l) and ( fT4l) . the two- 
parton distribution amplitudes can be expanded in a series of partial waves with different 
conformal spins. However, in analogy to that discussed in Refs. [H [9], by means of QCD 
EOMs, g^^\ g^^\ h^^\ and /ly^^ can be expressed in terms of twist-two two-parton and twist- 
three three-parton distribution amplitudes. The detailed results are given in Sees. 13.31 and 
13.41 It should be stressed that EOMs only relate terms with the same conformal spin so that 
the relations between LCDAs can be satisfied order by order in the conformal expansion [2j. 

3.2 Three-parton distribution amplitudes of twist-three 

The three-parton chiral-even distribution amplitudes of twist-3 are defined by 

{A{P, X)\qi{-x)-fal59sG^^{vx)q2{x)\0) = -p„[p^e£^ - p^e*}^^]f^^A{v, -px) + ... , (35) 

{A{P, X)\qi{-x)-fagsG^^{vx)q2{x)\0) = ipa[p^e*l^^ - P,y^±^^]f3A^{v, -px) + ... , (36) 
{A{P, \)\qi{-x)(Tai3l5gsGf,u{vx)q2{x)\0) 

= ~^^(^ [^"^/^^'i ~ P^f^aL - PaP^gp^, + PpPugi^]ftAfnAT{v, -px) + . . . , (37) 

where the ellipses stand for terms of twist higher than three, the following shorthand nota- 
tions are used: 

A{v,-px) = y"pae*^'"("«2-"'i+""«U(ag,ag„a3), (38) 
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etc., and the integration measure is defined as 



/T)a= / daq-^ / dag^ / dag (5(1 — 
Jo Jo Jo 



(39) 



Here a^j, ctgj, and are the respective momentum fractions carried by qi, q2 quarks and 
gluon in the axial-vector meson. Due to G-parity, for a l^Pi [1^-Pi] state A is antisymmetric 
[symmetric] under the interchange Ug^ in the SU(3) limit, while V and T are sym- 

metric [antisymmetric] (cf. the case of the p meson in Ref. [9]). Taking into account the 
contributions up to terms of conformal spin 9/2 and considering the corrections which arise 
from the quark masses, the distribution amplitudes can be approximately written as 

A^-i + o"3^Pi^(7ag - 3) , 



A = 5040(agj — aq^)aqj^aq^a'^ + SGOaq^^aq^ag 



V = 360aq-^ag^ag 



for the l^Pi states, and 



A = SGOaq-i^ag^ag 



1 + uj3p^-{7ag - 3)J + 5040(ag^ - a^Ja^y^a^^a^crap^, 
1 + uji-p^^{7ag - 3) + 5040(ag, - agjag^a^ja^aa^p^, 

1 + uj^p^^{7ag - 3) + 5040(ag, - ag2)ag,ag2ala^p^, 

1 



V = 5040(a;g^ — aq^)aqj^aq^a'^ + SGOaq-^^aq^ag 



T = 5040(agj — aq^)aqj^aq^a'^ + SGOaq^^aq^ag 



\^P^ + a^p^-{7ag - 3) 



(40) 
(41) 
(42) 

(43) 

(44) 
(45) 



for the l^Pi states, where A's correspond to conformal spin 7/2, while cj's and cr's are 
parameters with conformal spin 9/2. As the SU(3)-symmetry (and G-parity) is restored, we 
have A's=cr's=0. The normalization constants and are thus defined in such a 

way that 

j T^a{aq^ - aq^) A{ag^,aq^,ag) = 1, j Da A{aq^, aq^, ag) = A^p^, 

jvaV{aq^,aq^,ag) = 1, Jvaioq^ ~ aq^)V{ag^,ag^,ag) = aYp^, 

j 1^aT(aq^,aq^,ag) = 1, J Va{aq^ - aq^)T{ag^,aq^,ag) = atp^, (46) 



for the l^Pi states, and 



VaA{agj^,ag2,ag) = 1 



j Va{ag^ - ag^)V{ag^,ag2,ag 
Jva{ag,-aq,)T{ag,,ag,,a, 



1, J'DaV{aq^,aq2,ag) = XYp^, 
1, J Va T{aq-^,ag2, ag) = Xip^, 



(47) 



11 



for the l^Pi states. 



3.3 Relations among chiral-even LCDAs 

^\\{u,^) can be expanded in a series of Gegenbauer polynomials [HIS]: 



1=1 



(48) 



where /z is the normahzation scale and the multiphcatively renormalizable coefficients (or 
called Gegenbauer moments) are: 

= o.,^^^!w,^,^ o^ [' dxCf ^'i2x - l)^iix,^^). (49) 



3(/ + l)(/ + 2; 



i+i .\ 

+ 4$:^ , (51) 



In the limit of m^^ = rrig^, only terms with even (odd) / survive due to G-parity invariance 
for the l^Pi (1^-Pi) mesons. In the expansion of $|j^(M, /i) in Eq. (HHll . the conformal invari- 
ance of the light-cone QCD exhibits that partial waves with different conformal spin cannot 
mix under renormalization to leading-order accuracy. As a consequence, the Gegenbauer 
moments aj' renormalize multiphcatively: 

«" (/i) = aj' (/io) I ) ' (50) 

where b = (llA^^c — 2nj)/3 and the one-loop anomalous dimensions are 

with Cf = {Nl - 1)/{2N^). 

Applying the QCD equations of motion, as the case for the vector mesons in Refs. pfTO], 
one can obtain some useful nonlocal operator identities (see Appendix [B]) such that the 
two-parton distribution amplitudes g^f^ and g^^ can be represented in terms of and 
three-parton distribution amplitudes. Setting y = —x and adapting the formulas derived in 
Ref. P] for vector mesons to the present case, we fincifl 

[ due-'^P^g^f\u) = ! dt I rfMe-'*«P"<l)||(M) 
io Jo Jo 

-C^A^P^f I ^^^'^ / dvV{v,-tpx) - Ci^^ipxf [ dtt^ I dvvA{v,-tpx) 
Jo J-1 ' Jo J-1 

-hpxf [ dte [ due~''^P'^g^l\u) -i6+{px) [ dtt [ rfw e-**«P"$x(M), (52) 

^ .In .In .In .In 



^Here and below we use the notations close to that given in Ref. [9]. 
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and 



1 pi 

dtt / due-**«P^^^V) 
Jo 



— iC^j^{px) / dtt"^ dvA{v, —tpx) — i(^^^{px) / dtt'^ dvvV{v, 
Jo J-i ' Jo J-i 



-1 pi 

+ 'L I dtt due~''^f"'^±(u) 



'0 ^0 

where ^ = 2m — 1, and we introduced the abbreviation 



pV,A 

Ja ruA ' "'^'^ fArriA' 



Ia "^92 ± ^qi 



V,A _ JsA 



Solving Eqs. ( l52l) and ( l53l) . we obtain the sohitions [9j for g^l\u): 



U i 

g^l\u) = u jdv-'^{v)+u jdv-'^{v). 



and for g^^\u): 



dv - "^(v) + dv- ^(f 



- (5_<l> I (u) 



u u 



+ CIa [ ^"gi / ^"92 1 _^ (-r— + ) 

J J J- Q^gi '-^92 \ (^Olqi (l(y.q2 J 





M U 



.A d /■, ^(a) 







<Jl Clq2 



where 



+ 2C3% 







u u 



da„, I da, 



qi ^q2 



a 







1 ^^2 



doiq-^ dcX 



Oiq-^ (y.q2 



da 



A{a) 



V(«). 



92 , 
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Neglecting the three-parton distribution amplitudes containing gluons and terms propor- 
tional to light quark masses, g^^^ and g^^^ are thus related to the twist-2 ones by Wandzura- 
Wilczek-type relations: 

"1 



gy [u) ~ 



1 /■ 1 

dv -$||(f) + dv -$||(v) 



g^l^^^ {u) ^ 2u I dv-^,{v) + 2ul dv-^,{v) 



(58) 
(59) 



3.4 Relations among chiral-odd LCDAs 

The leading-twist LCDAs $^(u,/i) can be expanded as [H [2] 



$^(u,/i) = 6u(l - u) 



4 + E«^^('")^^'(2n-l) 



(=1 



(60) 



where the multiplicatively renormalizable Gegenbauer moments, in analogy to Eq. (149|) . read 

2(2/ + 3) 



dxCf/^(2x- l)$i(a;,/x), 



3(/ + !)(/ + 2) ,0 
which satisfy 

(/^) = {fW'"^) (/^o) 
with the one-loop anomalous dimensions being 



Qs(/io) 
a^(/i) 



(61) 



(62) 



(63) 



In the limit of ' with even (odd) / vanish due to G-parity invariance for the 

l^Pi (l^Pi) mesons. 

In analogy to Eqs. and (15^ . using the formulas developed in Ref. [H] for vector 
mesons to the present case, we obtain the integral equations: 



1 pi 



ipx / due-'^P''^h\ \u)-2 / due-'^'P'' {h\ '{u)-^±{u) 



= CsAiP^y I dvvT{v, -px) + {pxf [ e-'^P^^h^iu) + iS+px I ciw e-'^f^$|| (u), (64) 
J-i Jo Jo 

[ due-'^P''h^[\u) = -iC^^^px [ tdt [ dvr{v,-tpx)+ [ dt [ due-'^'P''h\'\u) 
Jo ' Jo J-i Jo Jo 



-6- f dt f due-'^"P''^\\[u] 
Jo Jo 



(65) 
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where we have introduced the notations 

= ^, Cs'yA = jr^ ■ (66) 

The solutions of the integral equations |9j given in Eqs. (l64ll and (165!) can be obtained to 



be 



u 1 

(1 + 6^)h\l'\u) =u jdv^ Q{v) + u jdv^ Q{v) , (67) 

u 



and 



h\f\u) = ijdv^&{v)-jdv^e{v) -6-^ii{u) 



u u 

^^SAI- /^"'72i ^(a)> (68) 

du J J 1- ttg^ - aq^ 





where 



en = 2<1.^(m) + ( <l>||(n) - ie<^>'|(")) -l^^'wiu) 



u u 



+ ^i^^ /^"- i-c,!-.,, + - ^''^ 



The two-parton twist-3 distribution amphtudes are thus related to the twist-2 ones approx- 
imately by Wandzura-Wilczek-type relations 

A'r "'(«) = 2 fa r + a /" dv^) . (71) 







3.5 Comparison of the vector mesons and axial-vector mesons 

In comparison with the vector mesons, we summarize in Tabled] the (anti-)symmetric prop- 
erties of LCDAs of axial-vector mesons in the SU(3) limit under the interchange of the 
momentum fractions of the quark and antiquark. On the other hand, comparing results for 
axial-vector mesons given in Eqs. d^HD-iEnD, dSS]), dMD, dSSD, and (IKTI - dRil) . and for 

vector mesons given in Ref . [9] , we can see an analogy between the two cases. One can relate 
the mathematical forms of two-parton LCDAs of twist-3 for vector mesons and axial- vector 
mesons in the following way: 
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vector mesons : axial-vector mesons : 

9± < > 9± 

(1 - « 

with the replacement 

+6± < — > -5=p 

+5± < — > -5=p 

A < — > V 

V < — > A 

r < — > r 

where the notations for LCDAs of vector mesons follow from Refs. [9], [in]. For vector mesons, 
additional factors (1 — 5+) and (1 — 5+) are considered [9] in the definitions of g^^^ and h\^^\ 

respectively, due to the normalizations g^f^{u)du = h^^^^du = 1. Nevertheless, we do 
not need to put such factors in the definitions of axial-vector mesons; see results shown in 
Eqs. (IMD, (I262D, ([263D, ([268]), ([269]), and ([270D. 

Table 1: The symmetric/antisymmetric properties of LCDAs under the interchange of the 
momentum fractions of the quark and antiquark in vector mesons and axial-vector mesons 
in the SU(3) limit, where the definitions of LCDAs of vector mesons follow from Refs. [9l[T0]. 



LCDA 


Vector meson 


Axial- vector meson (I'^-Pi) 


Axial- vector meson (1^-Pi) 


$11 


symmetric 


symmetric 


antisymmetric 


$± 


symmetric 


antisymmetric 


symmetric 


(a) 

91' 


symmetric 


symmetric 


antisymmetric 


{v) 

91' 


symmetric 


symmetric 


antisymmetric 


hf^ 


symmetric 


antisymmetric 


symmetric 


h\^^ or/i[|^) 


symmetric 


antisymmetric 


symmetric 


A 


symmetric 


antisymmetric 


symmetric 


V 


antisymmetric 


symmetric 


antisymmetric 


r 


antisymmetric 


symmetric 


antisymmetric 
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4 Determinations of leading-twist LCDAs 



4.1 Physical properties for axial- vector mesons 

In the quark model, the hght l^Pi states, referred to as 6i(1235), /;,i(1170), /ii(1380), and 
KiB, form the 1"^" nonets, whereas the hght l^Pi mesons, denoted as ai(1260), /i(1285), 
/i(1420), and Kia, form the l"*"^ nonets. /ii(1380) is not experimentaUy well-estabhshed [2^ 
and its quark content was suggested as ss in the QCD sum rule calculation |19]. It should be 
noted that the real physical states i^'i(1270) and i^i(1400) are the mixtures of l^Pi (Kia) 
and l^Pi (Kib) states; following the convention in Ref. [29], the relations can be written as 

|Ki(1270)) = \Kia) sinOK + l^^is) cos^k, 

|is:i(1400)) = \Kia) cos9k - \Kib) smOK- (72) 

In Ref. [29], two possible solutions with two-fold ambiguity |6'i<-| ^ 33° and 57° were obtained. 
A similar constraint 35° < \6k\ ^ 55° was found in Ref. [30j. Therefore, the favor values 
may lie in the range \6k\ — (45 ± 12)°. Just for simplicity, we will take 6^ = 45°, —45° 
as the reference points. Analogous to t] and rj', for I'^Pi states, /i(1285) and /i(1420) are 
mixed in terms of the pure octet fs and singlet /i due to SU(3) breaking effects, and can be 
parameterized as 

|/i(1285))) = |/i)cos^3p, + \fs)smesp^, |/i(1420)) = -|/i)sin^3p, + |/8)cos^3p, . (73) 
From the Gell-Mann-Okubo mass formula [3T1 [28] , it follows that 

cos' ^3p, = — ^ , (74) 

3 (^"^/l (1285) ~ "^/i(1420) ) 

where 

"^i^iA = "^Xi(i4oo) cos' 9k + m^,(i270) Sin' Ok ■ (75) 

Substituting into Eq. ([71]) with Ok = (45 ±12)°, we then obtaiifl^3pj = 38°+}6o which is con- 
sistent with the value of replacing m' by m throughout Eq. (17^ . The previous phenomeno- 
logical analyses suggested 6*3 p^ ~ 50° [32J. In the present paper, we will take 6*3 p^ = 38° or 
50° as the reference input for LCDA studies. 



Since 

4rn?^ — — "im?. 



tan e.p^ = '-^^^ > , (76) 



'KiA '"ai """'/l(1420) 

ns 

where m^g — (/ilTil/s) < with Ti. being the Hamiltonian, we find that 03p_^ > 0. In the present paper, we 
can extend the study in the traditional quark model to a field-theoretical consideration that each Fock state 
of /i [or fs] is proportional to {qq)"^g^{uu + dd + ss)/^/3 [or [qq)"^ {uu + dd — 2ss)/a/6], where there is a 
relative sign difference between the ss contents of /i and /s in our convention. 
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Similarly, for l^Pi states, /ii(1170) and /ii(1380) may be mixed in terms of the pure octet 
hs ad singlet hi, 

|/ii(1170))) = \hi)coseip^ + \hs)smeip^, |/ii(1380)) = sin^ip^ + j/ig) cos^ip, . (77) 

Again from the Gell-Mann-Okubo mass formula, we obtain 

3m? /117m — [4:711% —ml ) 

COS fiPi — -: r , [/is) 

•^y^hiiino) ~ "^/iiCiaso) ) 

where 

^K.B = ^Kr(1400) Sin^ 9k + m^i(l270) cos^ 9k . (79) 

We thus get 9ip^ ~ 10°^]^qo which coincides with the value of replacing by m throughout 
Eq. (!78|) . Note that I^kI > 50° is disfavored because of the constraint < cos^ 6*1 < 1. Note 
also that the QCD sum calculation suggested 6'ip^ ~ 45°, from which the predictive content 
of /;.i(1380) is predominated by the ss paii0 [19j. By comparing results with 9ip^ = 10° and 
45° in the phenomenological LCDA analysis, we will show the effects induced by different 
mixing angles. 



4.2 Axial- vector couplings of l^Pi mesons 

In this subsection, we calculate renormalizat ion-group (RG) improved QCD sum rules for 
the axial-vector couplings of l^Pi mesons. Because the QCD sum rule approach may not 
sufficiently determine the singlet-octet mixing angle, we thus calculate the QCD sum rules 
for pure l^Pi states. 



4.2.1 QCD sum rules for axial-vector couplings 

To evaluate the axial- vector couplings /ap^ for l^Pi mesons, we consider the two-point cor- 
relation function, 

= ^ J rf'a;e^^^(0|T(j^(a;)jt(0)|0) = ~Mq')g^, + U,{q')q^q, , (80) 

where the interpolating current = g27/^75?i satisfies 

(0|j^(0)|l=^Pi(P,A)) = -z/3p,m3p,eW. (81) 

In the massless quark limit, we have Hi = q^Il2 due to conservation of j^. Here we focus on 
Hi since Hi receives contributions only from axial- vector (^Pi) mesons, whereas 112 contains 
effects from pseudoscalar mesons. The above correlation function can be calculated from 

7 With the same reason as footnote [6l should be positive in sign. 
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the hadron and quark-gluon dynamical points of view, respectively. The lowest-lying l^Pi 
meson contribution can be approximated via the dispersion relation as 

4-%.if',.,i!^in£). (82) 

mip^ -q^ rr Jo s - 

where 11°^^, the QCD operator-product-expansion (OPE) result of Hi at the quark-gluon 
level, given in Ref. [19J up to dimension 6 and with 0{as) corrections, reads 



47r2 /i^ 



a. 



1 + — ) + 4 ( '^^) - m2{qiqi) - mi (5292) 



+ ^(^^«^(fe7MA"gi +g27MA''g2)5Zg7'^A"g) + 27ra,(giA"g2g2A"gi)Y (83) 

and Sq^^ is the threshold of the higher resonant states, such that the contributions of higher 
resonances are modeled by 



ds- 



vr I.^Pi s — q"^ 

We further apply the Borel (inverse-Laplace) transformation to both sides of Eq. fl82 

-1 n 



^4) 



B[/(g^)] = Jim (-g^) 



2\n+l 



(85) 



to improve the convergence of the OPE series and further suppress the contributions from 
higher resonances. Moreover, we adopt the vacuum saturation approximation for describing 
the four-quark condensates in the present work, i.e.. 



(0|gr*A'^ggr.A''g|0) 



1 



Tr(r,r,)Tr(A"A")(gg)^ 



(86) 



and neglect the possible corrections due to their anomalous dimensions. Finally, we arrive 
at the /apj sum rules, given by 



„9 9 -m? /M^ 



4^ ]s ds e-^/'^' (1 + ^) _ 1 {^G') + nig, {q,q,) + m,, {q,q,) 



M2 



((9292)^ + (gigi)^) + ^^q"' {q2q2){qiqi) 



(87) 
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4.2.2 Results for l^Pi mesons 



We start with the analysis of the /sp^ sum rules. To examine the quality of the sum rules, we 
also give the mass sum rule results. The mass sum rule for the l^Pi lowest-lying resonance 
can be obtained by taking the logarithm of both sides of Eq. ( IHTl) and then applying the 
differential operator M^d/dM"^ to them, where Sg^^ is determined by the maximum stability 
of the sum rules. Substituting the obtained Sq^^ and masses into Eq. (IHTI) . one arrives at 
the sum rules for the decay constants /a . In the numerical analysis, we use the parameters 
which are given in Appendix O and choose the Borel window 0.8 GeV^ < < 1.3 GeV^, 
where the contribution originating from higher resonances (and the continuum) is less than 
44% and the highest OPE term at the quark level is no more than 14%. Note that in the 
sum rules, the contributions from higher resonances are modeled by 



- / dse-^/*^'lmn^ 

^ J so 



OPE, X 



In Fig. [H the masses and decay constants are plotted as functions of the Borel mass squared 
M^. The results are summarized in Table [21 Introducing the decay constants /j^(i285) 

ffiiim 

(0|g7M75g|/i(1285)(P, A)) = -imj^^,285)fl^usA^^ ' (89) 
(0|g7M75g|/i(1420)(P, A)) = -^m^, (1420)/;, (1420) ' (90) 

we get 

//"(1285) = % ""^^ cosfep, + % ""^^ sinfep, = 173 ± 23 (167 ±22) MeV , (91) 

■"^ ' V 3 m/, (1285) V 6 m/, (1285) 

^..285^ = COS gap - ^^/? sin gap = -9 ± 13 (-59 ±18) MeV, (92) 

'^^(12^^) V^m/,(i285) ' V^m/,(i285) ' ^ ! ^ \ ! 

= -ll^ "^f' sinfep±% cos gap = -9±10 (-41 ±11) MeV , (93) 

^ ' V3m/,(i42o) V6m/,(i42o) 



//i(1420) 7f- Smgap, -=.- COS gap, 

V 3 m/, (1420) V6 m/,(i420) 

217 ±27 (-211 ±26) MeV, (94) 



corresponding to gap, = 38° (50°). In particular, /i(1285) and /i(1420) are predominated by 
their uu and ss contents, respectively. 

Using Eq. ( 1721) . the decay constants for -ft'i(i27o) and -ft'i(i4oo) read (with q = u or d) 



|g7^75s|iri(1270)(P, A)) = /i^i(i270) mxi(i270) eL^^ 



-i{fK^A^K^A sin dK + fK^s^K^Bf^Q^^" cos Ok) e^^\ (95) 
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Table 2: The sum rule results of masses, decay constants, and corresponding excited thresh- 
olds Sn^^ for l^Pi mesons. 



State 


Mass[GeV] 


Decay cosntant /spjMeV] 




ai(1260) 


1.23 ±0.06 


238 ± 10 


2.55 ±0.15 




1.28 ±0.06 


245 ± 13 


2.80 ±0.20 




1.29 ±0.05 


239 ± 13 


2.70 ±0.20 


KiA 


1.31 ±0.06 


250 ± 13 


2.90 ±0.20 



and 



(0|g7/.75s|^i(1400)(P, A)) = /i^i (1400) "^A-i (1400) eL^^ 



-i{fKiA^K^A cos Ok - fK.Bi^K^B^^o^^'' sin Ok) e 



.(A) 



(96) 



where fx^^^ '^k^a given in Tabled and ajj' (with f^^g = fk-^gi^ GeV) by definition 
in this paper) and mx^g can be found in Tables [3] and HI respectively. (See the detailed 
discussion in Sec. 14. 3[ ) Numerically, we thus obtain 



and 



/ki(1270) 
/i^i(1400) 

/Xi(1270) = 
/a'i(1400) = 



if taking Ok = 45°, 



197 ± 15 MeV , 
151 ± 12 MeV, 

-166 ± 11 MeV, 

179 ± 12 MeV , if taking Ok = -45° , 



(97) 
(98) 

(99) 
(100) 



where the correlation of the errors between the masses and decay constants is considered. It 
is interesting to note that, if one sets ajj''^^^ = 0, then the central value of /xi(i27o) becomes 
±182 GeV corresponding to 9k = ±45°, whereas that of /xi(i4oo) is 165 GeV independent 
of the sign of 9k', the local axial- vector current can couple only to Kia, but not to Kib in 
the SU(3) limit. 

4.3 Tensor couplings of l^Pi mesons 
The tensor couplings of l^Pi mesons are defined as 

(0|g2a^.gi|l'Pi(P, A)) = tf,\ e^,,/3ef,)P^ , (101) 
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Figure 1: masses and decay constants of the I'^Pi states as functions of the 
Borel mass squared, where the central values of input parameters given in 
Appendix [C] have been used. 
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{0\q2a''''l,q,\l'P,iP, A)) = -/i^, (ef.^P'^ - e^.^P^ • (102) 

In Ref. [18], we have calculated RG- improved QCD sum rules of the tensor couplings for 
l^Pi mesons, where we assumed that /ii(1170) is the same as 6i(1235), while /ii(1380) is 
made of ss. Here we instead examine the pure l^Pi octet states. The results are collected in 
Table [31 We have estimated the singlet-octet mixing angle by means of Gell-Mann-Okubo 
mass formula. Thus, the decay constants for the real states ^i(1170) and /ii(1380) can 
be obtained. See the results listed in Eqs. (11051) — (11081) . We have also updated the values 
for 6i(1235) and and Kib states due to the following two more reasons. First, we take into 
account the uncertainties of the condensates and quark mass which were not fully considered 
in Ref. [16|. Second, we re-examine the stability of the sum rules and choose to use the Borel 
window 0.8 GeV^ < < 1.3 GeV^ in the analysis. During this window, which covers the 
plateau region for the masses and decay constants versus M^, the contributions from excited 
states (including the continuum) and from the highest dimension term in OPE expansion 
are less than 43% and 5%, respectively, while for the previous choice with the upper bound 
1.5 GeV^ the contributions of excited states can reach 52%. Although the region for the 



present Borel window is a little lower than the previous study in Ref. [T8j, the present results 
should be more reliable. Note that it will become questionable for the sum rule results if 
one further reduces the lower bound of the Borel window. The reason is because in the 
lower Borel mass region the effects owing to the uncertainty of the highest dimensions term, 
for which we have assumed the vacuum saturation approximation, and radiative corrections 
become much more important and are out of control. 

Since the QCD sum rule approach may not sufficiently determine the singlet-octet mixing 
angle, we thus study the QCD sum rules for pure l^Pi states; in other words, we neglect the 
corrections arising from {hi\T-C\hs) due to {hi\T-C\hs) <^ {hi\T-C\hi), {hglHlhg), where TC is the 
Hamiltonian, in the studies of the mass and tensor coupling sum rules for hi,hs and Kib 
states. Employing the formulas given in Ref. pjj, we obtain the results shown in Table [31 
In complete analogy to the discussion for l^Pi states, we introduce 

{0\qcr,,q\h{im)iP, A)) = ^f,t;ilm W^fA)^'' , (103) 

(0|ga^.g|/ii(1380)(P, A)) = tfj^f^,,,,^ e,.^pe^,^P^ , (104) 

and then obtain 

fhSm = ^ cos Oip, + ^ sin ftp, = 116 ± 8 (128 ± 7) MeV , (105) 

/mii70) = ^ cos ftp, - ^ sin ftp, = 75 ± 8 (-36 ± 10) MeV , (106) 

/mi380) = ^'P^ + ^ ^'P^ = 58 ± 5 (-19 ± 7) MeV , (107) 
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/m1380) = -^'^^^^^1 ~ " ^ ^ ^ 

corresponding to 6ip^ = 10°(45°). As ftp^ = 45°, /ii(1380) is therefore dominated by the ss 
content, and /ii(1170) can be approximated by {uu + dd)/\/2. However, the ss content of 
/ii(1170) and uu content of /ii(1380) become significant for 6ip-^ = 10°. 
As for strange mesons, we have (with q = u,d) 

{0\qa,,s\K^{mO){P,\)) = z/i,(i270) W^rA)^"" 

= ^(/i,a^'^^^sin^;, + /^^^cos^^)e^,„^ef,)P'3 (109) 

and 

{0\qa,.s\K^{UOO){P,\)) = tf^^iuoo) W^fA)^'' 



^ifiA'""'" - fin e,uafsef,)P^ (110) 



where /^^^ and '^^'^ are given in Tables [3] and HJ respectively, and use of /^^^ = //^^^ = 
(250 ±20) MeV is made in the following by definition. Consequently, we obtain (at the scale 
/i = 1 GeV) 

/i(mo)= 145 ±15 MeV, (111) 
/ki(i4oo) = -124 ± 15 MeV , if taking Oj, = 45°, (112) 

and 

/x,(i270) = 124 ±14MeV, (113) 
/i(i4oo) = 145 ± 14 MeV , if taking Ok = -45°, (114) 

where ap '"^^^ gives about 8% corrections to the decay constants. If setting ad'^^'* = 0, all 
the magnitudes of the central values of the decay constants are equal to ~ 134 GeV. 



4.4 Gegenbauer moments of leading-twist LCDAs 

The Gegenbauer moments can be calculated from the standard QCD sum rule approach by 
adopting a relevant two-point correlation function as the starting point. It is interesting 
to note that in SU(3) symmetry limit the decay constant for a l^Pi state transiting to 
the vacuum via the local axial-vector current vanishes due to G-parity mismatch between 
the current and states. On the other hand, conventionally, the decay constants /sp^ for 
the local axial-vector currents coupling to the l^Pi states are chosen to be positive as in 
the present paper. Although the diagonal correlation functions may have good quahties 
for sum rule results of G-parity invariant parameters, it cannot determine the relative sign 
of the Gegenbauer moments. Moreover, the diagonal correlation functions cannot use to 
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Table 3: The sum rule results of masses, decay constants, and corresponding excited thresh- 
olds Sq^^ for l^Pi mesons. The values of f^p are given at the scale = 1 GeV. 



State 


Mass [GeV] 


Decay cosntant /f'^jMeV] 


So'' [GeV'] 


6i(1235) 


1.21 ±0.07 


180 ±8 


2.6 ±0.2 




1.23 ±0.07 


180 ± 12 


2.6 ±0.2 




1.37 ±0.07 


190 ± 10 


3.2 ±0.2 


KiB 


1.34 ±0.08 


190 ± 10 


3.1 ±0.2 



evaluate the G-parity violating parameters due to the mixing between ^Pi and ^Pi states. 
To determine not only the magnitudes but also the relative signs for Gegenbauer moments 
relevant to the leading-twist LCDAs of l^Pi states, we thus choose one of the interpolating 
currents in the two-point correlation function to be the local axial- vector current ^27At75Q'i, 
i.e., we consider the non-diagonal correlation functions here. Note that we define the sum 
rule to be diagonal here only if two interpolating currents in the correlation function are 
exactly the same, but non-diagonal otherwise (See footnote [9] for further discussions). With 
the same reason, in the following section we also adopt the local axial- vector current as one of 
the interpolating currents in QCD sum rule studies to determine the relevant parameters for 
twist-3 three-parton LCDAs of l^Pi states, whereas we use the local pseudo-tensor current 
as one of the interpolating currents since the resulting contributions arising from ^Pi states 
vanish in SU(3) limit. 

However, as we consider the non-diagonal correlation functions, G-parity breaking con- 
tributions of ^Pi and ^Pi states always mix. To obtain the G-parity violating Gegenbauer 
moments for LCDAs of ^Pi and ^Pi mesons, we will assume an additional reasonable con- 
straint on the parameters (see Eq. (11431) ). 

4.4.1 Gegenbauer moments of for I'^Pi mesons 

The LCDAs ^y^^ {u, fx) corresponding to the l^Pi states, which are denoted by the superscript 
^Pi, with the quark contents q2 and gi, are defined as 

{l'P^{P,X)\qM ^75g2(x)|0) = if,p,m-sp,{e<'h) rfxe^("P-^+^^'-^)$!^^(M, /i), (115) 

where u (or m = 1 — -u) is the momentum fraction carried by the quark qi (or antiquark ^2) 
and n is the renormalization scale of the LCDAs, and we have considered here and below 
that z (X y — X, z'^ = 0, {y — x)^ = 0, but z -/^ even for y ^ x. $||^^(m, /i) can be expanded 
in a series of Gegenbauer polynomials as given in Eq. fHHl) . To evaluate the Gegenbauer 
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moments of '^'y^S we take into account the following two-point correlation function 

n«(g) = z / d'xe^^^Oinn^lix) Ot(0)|0) = i-zqy[-I?iq') z, + I^\q'){zq)q,i (116) 



where 

O^(0) = g2(0)7M75gi(0) (117) 

and the relevant multiplicatively renormalizable operator, to leading logarithmic (LO) accu- 
racy, is 

I 

^3p,(a;) = Y,ci,,{tzd)'-^q2{x) {iz Oyqiix) , (118) 

j=0 

with D^=D^ - Df,= id +igsA%x)X''/2)^ - (d -zc/,A"(x)AV2)^ and Q,fc being the coeffi- 
cients of the Gegenbauer polynomials such that Cf^'^{x) = X]fc=o '^^fe^'^- ^sp^ satisfy 
the following relations: 

mi%mi'P^{P,X)) = -zMm3p,(6(^).-)(Pz)' ^^^ +^^^1^ + a^^f,) , (119) 

(0|O^(0)|l3Pi(P,A)) = -ifsp^msp/^^l (120) 

We are interested only in l'^\q^) since only states with quantum numbers of ^Pi contribute 
to lf\q^), whereas l2\(f') receives contributions from states with quantum numbers of 
pseudoscalar mesons and 3Pi. /P(g2) 

was already given in Ref. [33j. In the massless quark 
limit, one has if ^ = q^l'i'' ■ Nevertheless, the above relation is broken even for m^^ = ruq^ ^ 0. 
The OPE result of jj''', up to dimension 6 and with 0{pL^ corrections, is given by (for / > 1) 

,2 /•! 



if'''"' = -^q'ln^ [ duuuCf/\2u-l) In' - 



/^3/2/-|\ ^3/2 /-|\ 

- [m,Aq2q2)+m,M<li){-iy] + {-^G')[1 + (-1)'] 

+ Y^Cfll{l) 6(1 - l)[m,,{q^g,aGq^){~iy + m,,{q29sCrGq2)]L'''/^''^ 
167ras^3/2 



9g4 



cr^(l)(gi?i)feg2)[l + (-l)'] 



- [^f^'(l)]'[(^"i9i)'(-l)' + , (121) 

where the terms containing quark mass corrections and the results with odd I are new. For 
even / and in the massless quark limit, the above result for jf^^^ is consistent with q'lf^^^ 
given in Ref. 
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The quality of the sum rules obtained directly from jj'^^^^ is not good. See the discussion 
below. Another way to obtain the sum rules is to take into account the dispersion relation 
with one subtraction. This method was introduced in Refs. [Ml ES]- Consider if\q^) = 
~ '^^^^ {q"^) ■ ii\q'^) is finite in —q^ oo, and therefore we can write down 

^2 roo 



J«(g^) = /r(0) 



(0, 



(0) 



7o s(s-g2) 



ds 



[Pphys{s) - Im/| 



(I) ,pert I 



0] 



TT Jo s{s-q^) 



[p,.p,{s)-e{sl-s)lr^I?^^^^\s)], (122) 



where pphys and pisp^ are the total physical spectral density and the lowest-lying {l^Pi] 
spectral densit}|§, respectively, and where we have modeled the higher resonance states as 



Pphys{.s) = Pi3Pi(s) + e{s - s[l)ImJ| 



(0 ,pert/ 



(123) 



On the other hand, one should note that ^Pi states can still enter the sum results by the 
G-parity violating effect, which is due to m^^ — 7^ , 

{l'P,{P, \)MOh,lMm) = ^f^pA''^' rmp^ ef ) . (124) 

If one derives the Gegenbauer moment sum rules directly from Eq. (11211) . then the results 
are 



II 3p 2 2 -m? /M2 



27r3 



I % 



m2 /A/2 



3(/ + l)(/ + 2) 



2(2/ + 3) 

l-e->Yl + 4-1^)1 rduuuCf\2u-l) In^- 



Cf\l) ,«,^2 



12 



9M2 

for even I > 2 and 



22""^Cf/^(l) - ^'""^ r^3/.,,.,,_ , „ 



81M2 



(125) 



m\ /A/2 



3(/ + l)(/ + 2) 

2(2/ + 3) 



Cf'^{l) {mg^{q2q2) - mg^{qiqi)) 



+ JpCfll{l)e{l - l){mg,{q,gsaGq,)-mgM9saGq,))L-'"^'''^ 



+ 



(126) 



®Here we also need to consider l^Pi state which will contribute to the sum rules due to the G-parity 
violating effect. See Eq. p24|) . 
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for odd / (of the Kia meson), where 6{x) = 1 for x > or otherwise. However, if one 
chooses to use Eq. (11221) and divides the relation in Eq. (I122p by — before applying the 
Borel transform, then the Gegenbauer moment sum rules are 



3 Pi .2 



Pi 



27r3 



-m% IM^ 

e 



fi 



Pi 



3(/ + !)(/ + 2) 
2(2/ + 3) 



7(0 A 



■So 
M2 



1) In^ ^ 

u 



M2 



3/2 



1) 



12M2 ^Tl ' 



1 



+ 



'Cf\l){qiqi){m2) + 



9M4 

for even / > 2 and 



81M4 



[C^f^'(i)]'((gi9i)' + (g2g2)^) 



(127) 



(^l J 3 Pi 



e 



ll.'Pl^ll.'Pl ^^2 



Pi 



e 



3{l + l){l + 2) ^^,1 J, 
2(2/ + 3) 



M2 



mg,{q2q2) - mg^{qiqi)) 



^ Cf_/?(1) 0(/ - l){mg,{q,g,aGq,) - m,,(gi(?,aGgi))L-i^/(3^) 



81M4 



(128) 



for odd / (of the Ki^ meson), where we have substituted 



/f(0) 



^Pl r2 



3(/ + l)(/ + 2) 



2(2/ + 3) 27r3 



duuuCf'^{2u-l) In^ -, 

Til 



u 



which can be determined in the limit — > cxd, into the above two equations. 

Before we discuss the rehabihty of the sum rules about Eqs. (11251) . (I126P and (11271) . (I128p . 
one should note that, for even /, the contributions due to ^Pi states are ~ 0{rn?^ and thus 
negligible. Nevertheless, for odd /, the contributions for l^Pi and l^Pi states are of the same 
order of magnitude. The quality of the sum rules in Eqs. (11251) and (I126p is not good due to 
the following two reasons, (i) The OPE series converges very slowly. For instance, if taking 
a close look at Eq. (I125P with 1 = 2, the terms of dimension-4 and -6 in the OPE series are 
still comparable to the perturbative contribution even for choosing a quite large Borel mass 
~ 2 GeV^. (ii) The terms of dimension-4 and -6 have the opposite sign compared with 
the perturbative contribution, such that the contributions of higher resonance states that we 
modeled have also the opposite sign compared with the lowest-lying state. As a result, it is 
difficult to choose a reliable windows for sum rules given in Eqs. (11250 and (I126p . 
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Nevertheless, the sum rules, given Eqs. (11271) and (11281) . converge much more quickly. 
Meanwhile, in the OPE series in Eq. (11271) the contribution of the perturbative term have 
the same sign as the the terms of dimension-4 and -6. Consequently, we can find a suitable 
Borel window, where the contributions originating from higher resonances and the highest 
OPE terms are well under control. For the time being, we will hence focus on Eqs. (11271) 
and (I128p . The numerical results are given in Sec. 14.4.31 Considering Eq. (11281) with / = 1 
for the KiA (and Kib) mesons, the result approximately reads 



+ Cli do 



1 



-m2 /M2 



M2 



- ^ (m^gs^Gq^) - m,,{-sg,aGs))L~'^'^^') - ^ {{ssf - {^2?)^ (129) 
with m = {rriKiA + '"^i^is)/2- 

4.4.2 Gegenbauer moments of for l^Pi mesons 

To calculate the Gegenbauer moments of for l^Pi states, we consider the following 
correlation function 



1 / 

g^,u — -{qf^Zu + quz^) 

_ qz 



{qzY^'Uq') + 



(130) 



where 



a(0)=gi(0)7;.75g2(0). 



(131) 



and, to leading logarithmic accuracy, the relevant multiplicatively renormalizable operator 
is 

= Cn,,(^^5)'-^g2(x)a^„75^" {^z B)^gi(x) , (132) 



i=o 



with z being the light-like vector as defined previously, fi^''''' and O^, satisfy the relation: 
J](0|Oi(')(0)|l3Pi(P,A))(l^Pi(P,A)|O.(0)|0) 



ifipjsp^msp^ 



3(/ + !)(/ + 2) 
2(2/ + 3) 



Pz 



{P-zr'a\ 



+ 



Using the dispersion relation for Ti, a, ' ^ can be represented in the form 



3(/ + 1)(/ + 2) ftpjsp.m^p^ai 



2(2/ + 3) 



ml 



Pi 



1 To 

TT Jo 



'Pi 



ds 



OPE/ 



s — q'^ 



(133) 



(134) 
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where 



-i,c'i'\i)mq2) + {Mi){-ir') 



da[mq^a + mq^{a — l)\C^^'^{2a — 1) 



+ 



(1) + 2C;Z^(1)^(/-1) 



5/2/ 



3 
27r2 

3?6 



1+1 



2ocjii{mi - 2) + c^!:t{mi - 1)) (^G^)[(g2g2) + {q^K-iyn 

J TT 



.5/2, 



(135) 



Since ^Pi states can have smaU axial-vector couphng constants due the G-parity violating 
(SU(3)-breaking) effect, the RG-improved sum rules for Gegenbauer moments a/"' thus 
read 



±,3Pi r r± -mlp /M^ ^.ip^ ||.lp, ^ -m? /Af2 



3([+l)([+2)^_^x/, 



./ 3 ^2(l_g-4'^-VM^^ 



2(2/ + 3) 



{q29s(^ ■ Gq2) + {qigsO- ■ Gqi){-^) j-2/m 
M2 



TT 



2 r 



^Cf5(l)e(/ - 2) + ^cK(l)e(/ - 1)) 



x(^G^)((5-2?2> + (?-i?i)(-l)'+^)L^/^ , 

TT 



(136) 



where for odd / the corrections coming from the ^Pi states are of order and can be 
neglected, whereas for even / the I'^Pi and l^Pi states give contributions of the same order. 
With / = 0, 2, we can therefore obtain the following approximation for the K^a and Kib 
mesons: 



Oo ' + al 



X 



1 



'mKiAfx^A^KiA 'I^KiaIkia^Kia 



3 



M'\l - e-^°/^')(m,, - m,)L-^/^ - ((52^2) - (^s))L^/^ 



87r2 

_ {q2gsO' • Gq2) - {sgscr • Gs) 
3M2 



^"'^''^ + ^,{^G^)i{m2) - {ss))M, (137) 
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^2 ' + ^2 



with So ~ Sq 



4.4.3 Results 

In the numerical analysis, we shall adopt parameters which are collected in Appendix [Cl and 
Tables |2] and [31 It should be noted that for the moment sum rules the actual expansion 
parameter is M^/Z in the large I limit [?]. As a result, for a|' and af- with a larger / and fixed 
M^, the OPE series are convergent more slowly or even divergent as compared with the sum 
rules for the masses or decay constants. 

Consider af first. We adopt the sum rules given in Eqs. fll27p and fll28p . As discussed 
after Eq. fll28p . these two sum rules can give much more reliable results than those in 
Eqs. ([1251) and ([126]). For 

Qj2 1 expected, we find the higher Borel window to be 
2.0 GeV^ < < 3.0 GeV^, where the contributions originating from higher resonance states 
lie between 15% and 32%, and moreover the correction arising from the highest dimension 
term in OPE series is between 17% and 6%. For al'-^'^+a)'^^^a5'-^'^/|.^^//|:^^, the correction 
from highest dimension is quite small, but it is hard to estimate the contributions from higher 
resonances. Fortunately, the result of this sum rule is quite stable and we will choose to use 
the Borel window 1.5 GeV^ < Af^ < 2.5 GeV^, which is in between the cases of the sum 

II 3p 

rules for decay constants and 

The sum rules for ' ^ are the typical cases about the non-diagonal sum rules, for which 
the main contributions in OPE may come from the quark and quark-gluon condensates. It 
has been argued in Refs. [36], [371 [38] that these sum rules may suffer from contributions 
of higher resonances, so that the value about the lowest-lying meson may be overestimate. 
On the other hand, note that usually the radiative corrections are constructive at the 10% 
level for each OPE term. Since we have neglected the radiative corrections, it means that 
the real value for the lowest-lying meson may be underestimate. The above two corrections 
may partially cancel each other. Equivalently, these two corrections can be lumped into the 
uncertainties of the condensates and quark masses. In summary, the non-diagonal sum rules 
may suffer from above two corrections and the net effects on results may be less than 10%. 
Essentially, the above estimate is suitable for all twist-3 and twist-3 parameters studied in the 
present paper. However, since we do not do the qualitative estimates about the parameters, 
we thus do not include these possible errors in Tables [H [SI and [B 

The best stability of the af' sum rules is reached within the Borel window 1.5 GeV^ < 
< 2.5 GeV^, the same as the case of af'^^'^, where the correction from the highest 
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dimension term at the quark-gluon level is quite small. For G-parity violating Gegenbauer 
moments, Oq + C2 ■ a^'^^'^ and a^'^^^ + C2 ■ a^''^^-^aQ'^^^ with 

we find the suitable Borel windows to be 1.3 GeV^ < < 2.0 GeV^ and 2.0 GeV^ < < 
3.0 GeV^, respectively. Besides, the contributions originating from higher resonances (about 
15% ~ 32% and 15% ~ 32% for the former and latter ones, respectively) and the highest 
OPE terms (zero and 5% ~ 2% for the former and latter ones, respectively) are well under 
control. 

For G-parity invariant parameters involving quark mass corrections, because the results 
may be sensitive to the nonperturbative parameters for which we have considered larger 
uncertainties of the parameters here (see Appendix [C]), we therefore re-examine the errors 
of l^Pi states as well. The formulas for l^Pi states were given in Ref. [18]. The results for 
G-parity invariant Gegenbauer moments for l^Pi and l^Pi states are summarized in Table |H 
Note that the results for hi and hg (the l^Pi states) are new, where use of 
has been made in the numerical analysis and a^'"'^ = 0.0 ±0.1 [39j. To exhibit the quality of 
the sum rules for the l^Pi results, the Gegenbauer moments versus the Borel mass squared 
are plotted in Figs. [2] and [31 

For G-parity violating parameters, the Gegenbauer moments for Kia and Kib mix as 
given in Eqs. (1129^ . (I137P and (I138p . On the other hand, the G-parity violating quantities 
were not considered correctly in Ref. [18] owing to the mentioned mixtures and should read 



±,KiB , ±,KiA ±,KiA^^'KiAy 5 1 ^2 lj^j2 J ( 3 /- W 

- ^ {ms{q29sCTGq2) - m,,{sgsaGs))L-''/^''^ - ^(/^0'a^''*e-<•/^^' |, (139) 

and 

^KiB fkiB fKiB 18 '^KiB fkiB f^lB 

3 



87r 



_ ^^ {q29sa-Gq2)^-{sg.a.Gs) ^_,,,, _ |^ (^G^) ((,2,2) - {ss))L^^^] . (140) 

We use the updated values a^^*(l GeV) = 0.04±0.03 [IQ] and = (0.185±0.010) GeV [50] 
in the numerical analysis. Eqs. (11291) . fll37p . (11381) . (11391) . and (I140p do not offer sufficient 
relations to obtain explicit solutions for G-parity violating parameters but can give the 
following relations (at the scale ^ = IGeV): 
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ai'^'^ + (0.59 ± 0.15) al'^'^ = 0.17 ± 0.11 , 
_ (^37 ^ 0.56) a^'^'^ = 0.02 ± 0.08 , 
^±,KiA _ (0.01 ± 0.16) aj"^^^ = 0.02 ± 0.18 , 

^IK^B _ (^0.08 ± 0.07)ao '^1^ = 0.01 ± 0.09 , (141) 

where we have substituted the numerical values of mesons' masses and decay constants given 
in Tables |2] and [3l If the above G-parity violating parameters are expected to be small as the 
results for K and K* (for instance, see results in Ref. [IQ]), it is preferred that a^'^^'* and 
^ll.^s'is j^g^yg ^Y^Q same positive sign. On the other hand, comparing Eq. f ll36p in this paper 
and Eq. (3.23) in Ref. [IB], we obtain the good approximation between G-parity invarint 
Gegenbauer moments (with odd /) of l^Pi and l^Pi states: 



ar'.a'r'^^^^, (142) 

which can be reconfirmed from the results in Table HI Therefore, we further assume that 
the G-parity violating Gegenbauer moments satisfy a similar relation with an enlarged un- 
certainty range: 



,KiA 



"^K^A^K^Afk^A 



1.0 ±0.3. (143) 



From the above estimate, it follows 

'^lA _ n nc -I- n nn ^11'^ 



a^'"i^ = 0.08 ±0.09, a{l''^'^ = 0.14 ±0.15, (144) 
^IK,A ^ 0.00 ± 0.26, ai'^'^ = 0.17 ± 0.22 , (145) 



^2 



'^1-^ = 0.02 ± 0.20, a|'^'^ = 0.02 ± 0.10 , (146) 



which are collected in Table HI 

Finally, four remarks are in order. First, we will simply take f^p^ = fsp^, which is 

independent of the scale, in the study since only the products of ftp^af' are relevant. 
Second, the sum rules obtained from the nondiagonal correlation functions in Eq. (11301) can 
also determine the sign of ftp^a]^' relative to /sp^. Third, for the present case, the RG 
effects are relatively small compared with the uncertainties of input parameters. Fourth, 
neglecting the small isospin violation but considering the SU(3)-breaking correction, af,a| 
and are nonzero only for Kia- Note that here we adopt the convention that qi = s 
for KiA and Kip- For Kia (-^ib) containing an s quark, we have the following replacements 
'^15 '^0 2 ~^ ~(A^ ~'^o 2 ('^0 2? '^i" ~^ 25 "^i") G-parity violating parameters. 
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Table 4: Gegenbauer moments of and $|| for l^Pi and l^Pi mesons, where uses of 
fi~Pj^ — fspi and fip^ — fi~p^{i GeV) have been made. For l^Pi states, the results for hi and 
ha are new, for ajf and af'^''' are corrected, and for the rest are updated. 





II, m (1260) 

Z 


II f'^i 

Z 


Z 


„II,^1A 
tin 
Z 




1 GeV 
2.2 GeV 


-0.02 ± 0.02 
-0.01 ±0.01 


-0.04 ± 0.03 
-0.03 ± 0.02 


-0.07 ±0.04 
-0.05 ±0.03 


-0.05 ± 0.03 
-0.04 ± 0.02 


0.00 ± 0.26 
0.00 ± 0.22 


M 


±,oi(1260) 

i 




1 


1 


«0 


"2 


1 GeV 
2.2 GeV 


-1.04 ±0.34 
-0.83 ± 0.27 


-1.06 ±0.36 
-0.84 ±0.29 


-1.11 ±0.31 
-0.90 ±0.25 


-1.08 ±0.48 
-0.88 ± 0.39 


0.08 ± 0.09 
0.07 ± 0.08 


0.02 ± 0.20 
0.01 ±0.15 




II, 6i (1235) 

til 
i 


11,/^:^^ 


i 


i 


„ll,^ifl 
"o 




1 GeV 
2.2 GeV 


-1.95 ±0.35 
-1.61 ±0.29 


-2.00 ±0.35 
-1.65 ±0.29 


-1.95 ±0.35 
-1.61 ±0.29 


-1.95 ±0.45 
-1.57 ±0.37 


0.14 ±0.15 
0.14 ±0.15 


0.02 ±0.10 
0.01 ± 0.07 


M 


±,6i(1235) 
"2 


^2 


^2 






1 GeV 
2.2 GeV 


0.03 ±0.19 

0.02 ±0.15 


0.18 ±0.22 
0.14 ±0.17 


0.14 ±0.22 
0.11 ±0.17 


-0.02 ± 0.22 

-0.02 ±0.17 


0.17 ±0.22 
0.14 ±0.18 



5 Determinations of three-parton LCDAs of twist-three 

In this section, using the QCD sum rule approach, we estimate the relevant parameters 
involving SU(3)-breaking effects in determinations of the three-parton LCDAs of twist-three. 
The quark masses (SU(3)-brcaking effects), can give contributions not only to G-parity 
invariant conformal moments but also G-parity violating ones of LCDAs. In the calculation, 
the OPE series in the QCD sum rules are evaluated up to dimension-seven, but up to 
dimension-four for terms proportional to the quark masses. 

The parameters can be calculated from the standard QCD sum rule approach by adopting 
a relevant two-point correlation function as the starting point. The choices of suitable 
interpolating currents may affect the qualities of the final estimates. Note again that in 
SU(3) limit the decay constant for a l^Pi state transiting to the vacuum via the local 
(pseudo-)tensor current vanishes due to the G-parity mismatch between the current and 
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Figure 2: G-parity invariant Gegenbauer moments of leading twist LCDAs for 
l^Pi states, corresponding to the scale fi = I GeV, as functions of the Borel 
mass squared. The solid curves and bands correspond to the central values 
and uncertainties of the input parameters, respectively. 
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Figure 3: G-parity violating Gegenbauer moments of leading twist LCDAs, 
corresponding to the scale = 1 GeV, as functions of the Borel mass squared, 

where Ci = fj^s/fkA ^2 = {mK,s/mK,^){f^^^fK,s/fk^Afi<iA)- The solid 
curves and bands correspond to the central values and uncertainties of the 
input parameters, respectively. 



states, whereas the axial-vector decay constant vanishes for a l^Pi meson. Motivated from 
the above properties and to determine not only the magnitudes but also the relative signs for 
parameters relevant to three-parton LCDAs of twist-3, we thus choose one of the interpolating 
currents in the two-point correlation functions to be the local axial-vector (or pseudo-tensor) 
current to calculate the parameters for the l^Pi (or l^Pi) state, i.e., we consider the non- 
diagonal correlation functions here. We did not consider the diagonal correlation functions 
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in the beginning since it cannot determine the sign □ of the parameters and cannot be used 
to evaluate the G-parity violating parameters due to the mixes of the l^Pi and l^Pi states. 

Unlike results from diagonal correlation functions, where the perturbative contribution 
may dominate in the OPE series, our results show that the leading contributions are dom- 
inated by the term with the quark or gluon condensate in the most cases, and we may 
need to take into account a higher Borel window, so that the contribution originating from 
the highest dimension term in the OPE expansion can be well under control. However, a 
problem occurs in the /^sp^ study for which the quark condensate contributions vanish in 
0{as) after adding all the diagrams, and therefore the OPE result become highly unreliable 
(see the details in Sec. I5.1.3p . We therefore resort to the diagonal sum rule for pursuing 
this parameter. However the lowest lying pseudoscalar meson contributes to the diagonal 
sum rule. Fortunately, although the pseudoscalar meson contribution is involved in the di- 
agonal sum rule, such effects can be subtracted by using the non-diagonal sum rule for the 
pseudoscalar meson. These two sum rules have been studied in Ref. pOj for calculating the 
twist-3 parameter f^K relevant to the twist-3 three-parton LCDA of the kaon. Instead, we 
use these two sum rules to extract the value of f^sp^- 

We study the parameters relevant to twist-3 three-parton LCDAs for l^Pi and l^Pi 
states in Sees. 15.11 and [5^ respectively. One should note that, as calculating the sum rules 
for G-parity symmetric parameters of twist-3 three-parton LCDAs of l^Pi (1^-Pi) mesons, 
the G-parity violating parameters for l^Pi (1^-Pi) states, relatively suppressed by O^m"^), 
contribute to them but can be negligible. However, the sum rules for G-parity violating 
parameters for l^Pi (1^-Pi) mesons receive the contributions arising from G-parity symmetric 
parameters about l^Pi (l^Pi) states, which are of the same order of magnitude because the 
axial-vector decay constants (pseudo-tensor decay constants) for l^Pi (1^-Pi) states do not 
vanish due to the quark mass corrections. The numerical results are given in Sec. 15.31 To 
LO approximation and including the quark mass corrections, the complete RG evolutions of 
the parameters are given in Appendix O 

5.1 Axial- vector mesons with quantum number l^Pi 
5.1.1 /g^apj, t^apj and aYp^ 

The coupling constants /^sp^, ujYp^ and aYp^ can be obtained by considering the correlation 
functions, 

1 1 d'xe^'^^ {0\t{j!'X{x), gi(0)7.75g2(0)}|0) = -rrp^iq') {qzYgi^ + ■ • • , (147) 

^ It is a little different for the definition of the diagonal correlation function here and in Ref. ^U\. In 
Ref. [2^ and related studies, the authors choose the non-local light-ray operators in calculating the correlation 
functions. According to their definition, although some parameters' signs can be determined in the diagonal 
sum rules which are actually non-diagonal in our definition, the relative signs of f^A^fs^A fi'A ^^il^ 
cannot be established. 
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i J d'xe"^^ {0\t{j'Z{x), ?-i(0)7.75?2(0)}|0) = -T^^^^iq') {qzfgi^ + ■■■, (148) 

i j d'xe^"^ {0\t[jI;1{x), gi(0)7.75g2(0)}|0) = -T^^M') ^^I'^fdtu + • • • , (149) 
where the currents are defined as 

■/li (0) = ^"^^?2(0)7a^.G^M(0)gi(0) , (150) 

4:^(0) = z^z^U^)lags[iDzGp^mqM , (151) 

4:^(0) = z'^z^HOhagstz'iGfs^iO) Ds - Ds G^^(0)]gi(0) , (152) 
which can couple to the l^Pi states as 

{0\j!:lmi'P,{P, A)) = -zfl^p^iP-zfe^ll + 0{-z,) , (153) 

{0\4'Xmi'P,{P, A)) = -zf^sp^ia^KP-zfe^ll + 0{z,) , (154) 

{0\j!:lm^'P,{P, A)) = -zflsp^aYp^iPzMl + 0{z,) , (155) 

with 0{z^) including the twist-4 correction and the average gluon momentum fraction {a^) 
satisfying 

= ^ + ^< ■ (156) 

Here and below the ellipses denote terms irrelevant to the present studies. It should be note 
that T^p^, Tip"", and T^p^ receive no contributions from pseudoscalar states. 

Assuming the quark-hadron duality, we can obtain the approximate expressions 



1 y 1 f ImTsp [s) 

— ^mzpjzpj^zp = - ds , (157) 

mip^ -q^ ' n Jo s - q^ 

—2 3^^PiM/3mK) = - / ^^-^ , (158) 

"hp^ ~ q Jo * ~ y 



and 



where 



1 V V 1 r° Imr3^'^'°P^(5) 
— j"^3Pi/3Pi/3,3Pi(^3p^ = - / ds , (159) 



1447r3^ ^2 247rg2 ^ ^ ' ' 27q 
a 



367rg2 



5(mi(gigi) +m2{q2q2) ) +4(7711(52^2) +r"2(gigi) ) ( In — |- - — 



At^ 6 

(160) 
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y,ag, OPE 



3Pi 



r^'ln ^ + (0 + 0{as)){asG') - {mi) {^2) 



' ^ (mi{qiqi) + m2{q2q2)^ + (mi{q2q2) + m2{qiqi)^ ( In 



-g2 11 



(161) 



and 



^ ^mi{qiqi) -m2(?2g2)j - (mi{q2q2) - m2{qiqi)^ (^\n^- - y 



+0 • {asG') + • (gigi)(g2?2) + (0 + C»(a,^)) (gsgs)') ■ (162) 

Note that T^p_^, T^^, and T^p^ receive contributions from ^Pi states because ^Pi states have 
small axial- vector coupling constants due to ruq^ — niq^ 7^ 0. Consequently, after performing 
the Borel transform, we obtain the sum rules: 



-ml /M^ V , -ml /M^ ||,liPi ,y 

e mspJspJlsp^ + e mipjip.a'^' flip^Xli 

5(rni{qiqi) + m2{q2q2)^ 



367r 



-4 (mi (5252) + m2 (qiqi)) + 7e + In ^ - Ei - ^ 



(163) 



/M2 



28 



1,0 /jvj c cV / 3 3 y \ — m? /M^ II l^PwV / 3 T/ 

e mspjapjgsp^l - + — o^spj +e ^1 mipjip^al hipA^^ip, 



1207r3 



1 1 ^ 

+ (mi{q2q2) + m2{qiqi)) ( "3" + + ~ 



So 
M2 



(164) 



e ^1 mspjspj^sp^aap^+e ^1 mipjip^a'^ /gip^ 



367r 



-(mi(gigi) -m2{q2q2) 



+ (mi (52^2) - m2(gigi)) + + In - Ei (^ - 



(165) 
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As for Eqs. 



and (11641) . the contributions, arising from G-parity breaking parameters 



relevant to twist-3 three-parton LCDAs of the l^Pi state, are relatively suppressed by 0{rn?^) 
and can be negligible. Nevertheless, in Eq. fll65p . the contribution originating from the l^Pi 
state is of the same order of magnitude as the G-parity breaking parameter relevant to the 
l^Pi state, and should be taken into account. Similar situations occur in Eqs. (11781) . (I179p . 
(fTSOD . (fT93D . (fTMD . and (fT95|) in the following subsections [5X1 and [5X1 

5.1.2 /3^3pj, A^^, and aip^ 

The constants /^ap^, A^^, and aip^ for l^Pi mesons can be defined as the following matrix 
elements 



(0| JJ^(0)|A(P, A)) = -/3^3p, {Pzfe% + 0{-z,) , 

(0| jJ;t(0)|/l(P, A)) = -/3%^ (^A^ + ^ai)i {P-zMl + 0{-z,) , 
where 0{z^) involves the twist-4 correction and the interpolating currents are 

= z"z^q2{^)iai,gs[iD-zGp^{Q)\qM , 

= ^''^^Q2{0hal59s[Gf,,{0)i Dz-tD zGp,m(lM ■ 
Thus, to evaluate f^sp^, A^sp^, and c^sp^, we consider the non-diagonal correlation functions. 



(166) 
(167) 

(168) 

(169) 
(170) 

(171) 



d'xe^^^ (0|t{ Jj^(a;), gi (0)7.7552(0) } |0) = T,^p^{q') {qzfg 
d'xe^^- (0|t{ Jf^(x), gi(0)7.75g2(0)}|0) = T.'^^.^iq') (qzfg 
d'xe^^- (0|t{ Jj^(x), gi(0)7.75g2(0)}|0) = T.^^f^iq') {qzfg 



+ 



+ ■ 



+ ■ 



(172) 
(173) 
(174) 



where we get 



T3 



as 2i._ 
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miiqiqi) + m2{q2q2) 



2 „ \2 



+ i0 + O{al))({qiqi){q2q2), {qiqi)\ {^2) 



(175) 
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5(rni{qiqi) - m2{q2q2)^ -4^mi(?2g2) - m2{qiqi)^ ( In 



2 



3 



+o-(«.G^)-|^((g-igi)^-(g-2?2)^), 



(176) 



yA,o-,OPE^^2 



247rg2 



miiqiqi) - m2(g2?2)) - 2(^mi(52?2) - "T,2(gigi)) f In 



+0 • (a.G^) + • (gigi) (?2?2) + (0 + (?-2g2)') ■ (177) 

From the above OPE results, it follows the sum rules: 

-m'i„ /M^ „ „4 -m? /M^ „ WI^Pt rA A 



14407r3 



«5 / ° _ _-s/M2 



727r 



367r 



mi{qiqi) +m2{q2q2)^ 



+ (mi(^2?2) + m2{qiqi)) Q - 7b - In + Ei - -^^^ 



(178) 



e '-^1 



({miY - {M2f^ + ^ - 5(^mi(^igi) - m2(?2?2)) 

2 



'27M2 



(179) 



and 



^mi(gigi) - m2{q2q2)^ 



3 3^ 
28< 



+2(^7711(52^2) - 7772(gigi)) - 1 - 7e - In + Ei|^ - 



(180) 



5.1.3 f^sp^, ujtp^ and asp^ 

To evaluate the couphng constants f^sp^ , oj^p^ and cr^j'^^ of l^Pi states, we first consider the 
following non-diagonal correlation functions, 

7 j ci'^a;e^^^(0|T{j^^(a;),gi(0)7^75g2(0)}|0) =T3i,^(g2)(g^)^^ + --- , (181) 
i j d^a;e^''-(0|r{4'^(x),gi(0)7^75g2(0)}|0)=T3-i;;^«(a(g^)% + -- - , (182) 
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1 1 d'xe''^^ (0|t{j3^'^(x) ,gi(0)7^75g2(0)}|0) = T^^^iq^) {qzfz, + ■ ■ ■ , (183) 

where 

^i''^(0) = g2(0)^%Aa„^75^7sG""(0)gl(0) , (184) 

j3'^(0) = g2(0)^^^Aa«;375^7.[^^^G'^"(0)]gi(0) , (185) 

^3'^(0) = g2(0)^%Aa„/375^7.^^'[G^"(0) Ds - Ds G'^"(0)]gi(0) . (186) 

To including the hght quark masses consistently to the 0{as), the currents needs to be 
replaced by the renormalized ones as 

^ = + ^±.(mi+ m2)z''tD ■ z{q2 ^.qi) , (187) 
47r ye 

j3.^ — . J3'^ = + ^l(m, + m2)z^zD ■ z)\q2 ^l,qi) , (188) 

^3'^ — ^3'^ = ^3'^ - ^^(^1 - m2)z'^i^D . -zf{q2 ^75gi) , (189) 



where 



11 

- = — h 7_B — ln47r . 

e e 



and d (dimension) = A + 2e. As adopted in all the calculations of this paper, the modified 
minimal substraction scheme (MS) is used to regularize the divergent integrals. From the 

3 _L 

technical point of view, the reason that we have to take into account the mixings of ' and 
twist-2 operators is because we need to remove the nonphysical ln(— g^//i^)/e terms in the 
calculation. Physically speaking, /^ap^ and f^zp^ojtp^ mix with fip^ijUq^+niq^), while f^^p^atp^ 
mixes with fzp^{mq^ — rriq^). For a massive quark, the above mixings were studied in Ref. 
[20] by using the light-ray-operator technique |41j. According to their results, (i) /^3p^ can 

mix with (m,^ +mq^) f^p^ and (m,^ +mq^) /sp^af' ^\ (ii) f^^p^ujtp^ with (m^^ +mq^) fsp^ and 

(mg, + 771,2 ) /3Pi«i' and (m,^ + niq^) f-ip^a^ ^\ and (iii) f^^p^atp^ with ( 

) /sPi 

and (m,^ — m^j) /spjaf' and (m,^ — m^j) /spjaf Neglecting the corrections arising 

II 3p 11 3p 

from a^' and , our results basically agree with those given in Ref. [20] except a 

3 _L f 

sign difference for Jg' . To LO approximation, the complete RG evolutions of the relevant 
parameters are collected in Appendix ICl 

Tg-^^ (g^), Tg^"", and T^'^ are relevant for the present consideration and the OPE results 
are given by 

Tt^''''{q') = (mi+m2)^gMn^[5-21n^] +(O + 0(a^))((gigi) + (g2g2)) 



/i 

1 /-q^ 547\a 



( In -J- - — ) ^ {{qigs(yGqi) + {q29s(^Gq2)^ 
— {a,G^){{q,q,) + {q,q,)), (190) 
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(191) 
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847\a 



96 / TT 



— (gi^f^aCgi) - {q2gs(^Gq2) 



)((gigi) - (g2g2))- 



(192) 



For the coefficient of the quark-gluon condensate in the above equations, the logarithmic 
factor arises from the contribution of the diagrams that contains a loop connecting a quark 
propagator, where the quark propagator (which does not belong to the part of the loop) emits 
a soft gluon into the condensate. Note that there is no infrared pole (IR) in the calculations 
of the present work due to the fact that the off-shell external momentum —q^ < regularize 
the IR singularity in QCD sum rule approach. (It is interesting to note that the infrared 
sensitive terms ~ ln(— g^/m^) may appear as considering the order up to m^. However, they 
can be absorbed into the condensate [42j.) 

The sum rules for f^sp^, wg^^ and cr^'^^ therefore read 
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and 



Imi - m2, 



57671^ 



5 /i^ 



+ 7^((9i?i) - (92^2)) 
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a. 



2167r 

vr 
9M2 



847 
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M2 



{qigsCrGqi) - {q2gsCrGq2) 



(195) 



Unfortunately, one can read from Eq. fllQOp or Eq. fll93p that for the f^sp^ sum rule the 
term of dimension-3 in OPE series vanishes in 0{as)i while the terms of dimensions-5 and 
-7 are comparable in magnitude but with the opposite signs. As a result, the OPE series 
does not show any convergent behavior. 

Therefore, to evaluate /^sp^, we further consider the following diagonal correlation func- 
tion, 

(196) 



rf^xe'"- {Q\T{ jf ^^(x) , Jf ^^(0) = -f\q') {qzf . 



However, at the hadronic level, the lowest-lying resonance for the above correlation function 
is the pseudoscalar meson (the 0~ state). To subtract the contribution arising from the 
lowest pseudoscalar meson, we study the following non-diagonal correlation function: 



where 



d'xe'^- {0\T\ J^f H^) , gi(0)75g2(0)||0) = -T'\q') {qzf 



z^z^^{^\q2a^pl^gsG^^qi\PS{P)) = 2ihps{Pz) 



(197) 



198) 



with "PS* = the lowest-lying pseudoscalar meson" . Only the pseudoscalar mesons contribute 
to the above non-diagonal correlation function. T-^{q^) and T' (g^) have been calculated in 
Ref. |20] for studying f^^K- 

Using the results given in Ref. |20j, the corresponding sum rules are 



89 as 
5184^ 



+ Y^("^gi(^i?i) +"^92(^2^2)) - 103^ jvf^ [mqAQi9s(rGqi) + mq^{q2gs(yGq2) 



729 M2 
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(199) 
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and 



,PS 
'0 



-^("^gi (^2^2) + mg, {Ml)) (^3 + + 1^ M2 - ^1 ( - M2 
^niq^ {q2gs(^Gq2) + nig^ {qigsCrGqi) 



6M2 



16 TTO^ X6 TTQJ 

+ ^j^iiiiqi)' + (9292)') + Y^(gigi)(g2g2) • (200) 

In calculating the f^sp^ sum rule, we substitute f^ps by using the expression given in 
Eq. f l200p . Although the sign of fj^sp^ compared with fsp^ cannot be determined from 
Eq. fll99p . there are two indications that the sign of /^sp^ should be negative. One is that at 
the large limit the result of Eq. (11931) implies the negative /^sp^ although the sum rule 
cannot offer its reliable magnitude. The other one is that only the negative /^sp^ can result 
in a physical value for (aj"), the average gluon momentum fraction in a l^Pi meson, which 
should satisfy 

< «) = ^ + 1^3^, < 1 . (201) 

Note that, as for cj;^^ sum rule (Eq. (11941) ) which is dominated by the quark-gluon condensate, 
because the contribution 0{as) of the dimension-7 term, ~ {c(sG'^){{qiqi) + {q2q2)), is absent, 
we thus do not further study a different correlation function. 



5.2 Axial- vector mesons with quantum number 1 Pi 
5.2.1 fl^p^, XYp^, and aYp^ 

The coupling constants f^ip^, XYp^, and aYp^ can be obtained through the following matrix 
element 

(0| JS(0)|1^Pi(P, A)) = -./3%, {P-zMi + 0{-z,) , (202) 
(0| JK(0)|1^Pi(P, A)) = -./3%^ XYp^P-zMl + 0{z,) , (203) 

(0| JK(0)|1^Pi(P, A)) = -^fl.p^ {^-XYp^ + ^^aYp}j {PzMl + 0{z,) , (204) 

where the currents have been defined in Eqs. (I150p . (I15ip . and (I152p . To evaluate f^\p^, 
X}Yp^i and uYp^i we consider the following "non-diagonal" correlation functions, 

j rf^xe^^^(0|T{jJ';:(x) ,gi(0)zV.A75g2(0)}|0) = TYp^{q^) [qzYg^, + ■ ■ ■ , (205) 
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d'xe^^^{Q\T{4'^^{x),q^iO)z'a,,^,q2{0)]\0) =Tr^^^^^^ (206) 

d'xe'^^0\T{4;l{x),qm^'a,xl,q2{0)]\0) = TYp^^) {qz)'gi, + ■■■. (207) 

Concerning the light quark masse corrections to the correlation functions, we have to replace 
the currents by the renormalized ones: 

Jt, — ^ J!'^, = JZ + £t^(^i + ^^r^'^^i^D . -zm2<ya,l,qi) , (208) 

JZ JZ = JZ - ^^(^1 - m2)z'^titD ■ -z){q2^a,l,qi) , (209) 

JZ, — JZ, = -JZ - ?^T^(^i - m2)z'^^i^D . -zma^,^,q,) . (210) 



The above mixings lead to that f^ip^ mixes with f^p_^{mq-^+mq^), while fYip_^Xi'p_^ and fYip_^(T^_^ 
mix with f^p^{'mq^ — niq^). We did not find any explicit result in the literature that can 
be used to compare with the present calculations. However, in analogy to the discussion 
after Eq. (I189p . the relevant parameters can mix in addition with f^p_^{mq-^ ± 771^2)^1"' 
and f^p^ijnq^ ± m^Ja^' where the upper sign corresponds to the G-parity conserving 
parameter f^ip^ and the lower sign to the G-parity violating parameters for which f^ip^^ip^ 
does not mix with fi~p_^{mq-^ — m^Ja^' We have neglected the RG-corrections due to 
(mgj ± m^Ja^' and (m^^ ± m^Ja^' in the present calculations. Considering the mass 
corrections in the RG equations, the scale dependence of the parameters relevant to the 
twist-three three-parton LCDAs is summarized in Appendix [Cl 
The OPE results of TYp^{q^), TYp\{q^) and Ti^"(g^) are 
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respectively. Note that TYp^, Tip^ , and can receive contributions from ^Pi states because 
^Pi states have small pseudo-tensor coupling constants due to the unequal quark masses. 
Consequently, we obtain the QCD sum rules 
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In Eqs. (12141) . the contribution originating from G-parity breaking parameters a^'^ ^^"^sp^ 
relevant to LCDAs of the l^Pi state are relatively suppressed by 0{rn?^ and can thus be 
neglected. In Eqs. (12151) and (I216p . the contributions arising from the l^Pi state are of the 
same order of magnitude as the G-parity breaking parameters related to the l^Pi state, and 
should be taken into account in the numerical analysis. Analogously, one needs to take into 
account the corrections due to the l^Pi state in Eqs. (12321) . (I243p and (I244p in the following 
subsections, whereas such corrections can be negligible in Eqs. (I230p . (I23ip and (I242p . 

5.2.2 f^^p^, u^p^ and a^p^ 

The parameters f^ip^, ^ip^ and af^^ are defined through the matrix elements: 
(0| jf;(0)|l^Pi(P, A)) = -f,%,{P-zMl + 0{z 
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(217) 



(0|Ji^(0)|l^Pi(P,A)) = -/3%,K)(P^)^e^i + 0(^,), (218) 

{0\Jl:imi'Pi{P, A)) = -/3^.p,<(P^)^e^], + 0{-z,) , (219) 

where the interpolating currents have been given by Eqs. (11691) . fll7Up . and (11711) . 0{z^) 
contains the twist-3 and twist-4 corrections, and the average gluon momentum fraction (a^) 
satisfies 

«) = ^ + |< • (220) 
/^ip^, u^p^ and a^p^ can be therefore evaluated by considering the correlation functions, 

1 1 d'xe^'^^ {0\t{j!:^^{x) , gi(0)zV,,75g2(0)}|0) = T,\{q') {qzfg^^ + • • • , (221) 
1 1 d'xe^^- (0|t{ J|;;l(x) , gi(0)zV.,75g2(0)}|0) = Tf^^q') {qzYg^, + ■ ■ ■ , (222) 

I j d'xe^'^^ , 9i(0)^V.A75g2(0)}|0) = T,^p;(g2) {qz)Xu + " " " , (223) 

respectively. To consider the light quark masses consistently to the 0{as), the currents needs 
to be replaced by the renormalized ones as 

■J'd Jli = + + ^2)n^D ■ z){q2<y.,imi) , (224) 

4^ = + + ^^Yz'-iiD ■ zf{q2<y.,l,qi) , (225) 

^3,^ — ^3,'^ = Jli - - ^2)n^D ■ zf{q2<ya,im) . (226) 

From the above results we obtain that f^ip^ and f^ip^ujip^ mix with fi~p_^{mq^ + m^j), while 
fs^Pi'^^Pi iiiixss with f^p^iruq^ — nig^). As our results in the previous subsection, we did not 
find any literature that can be used to compare with the present calculations. Again, in 
analogy to the discussions after Eqs. (11891) and (I210p . the relevant parameters can mix in 
addition with f^p^{raq^ H-rngJa^' for /^ip^, and with f^p^{raq^ im^Ja^' and f^p^{mq^ ± 



niq^ )a. 



^' ^'^ for f^ip_ujip^ corresponding to the upper sign and f^ip^a^p^ to the lower sign. We 

neglect the RG-corrections due to (m^^ ± mq^)af' and (m^^ ± m^Ja^' in the present 
calculations. The RG evolutions, containing the quark mass corrections, for the relevant 
parameters are summarized in Appendix O 
We get 
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The resulting QCD sum rules read 
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5.2.3 /g-^ip^, Xip^, and a^p^ 

The coupling constants fj^ip^, Xip^, and a^p^ for l^Pi mesons are defined as the following 
matrix elements 

(0|J3^'^(0)|l^Pi(P,A)) = ./3^pjP^)^(eWz), (233) 
(0|J^^(0)|liPi(P, A)) = zf,\p^ Afp^(Pz)(6Wz) , (234) 

(0|J,^'^(0)|liPi(P,A)) =./3^P^ (jXlp^ + ^alp^jiP-zn^^'^-^), (235) 

where the interpolating currents have been given by Eqs. fll84p . (11851) . and (I186p . To 
calculate these three parameters, we consider the following correlation functions, 

J d'xe"^^0\T[4^{x), gi(0)zV^,75g2(0)}|0) = -T.'p^iq') {qzfz, + ■■■, (236) 
J d^xe"^^{0\T[jf^{x), gi(0)ra^,75g2(0)}|0) = -T,-^pf (g') iqz)% + ■■■ , (237) 
d^xe"^^{0\T[4^ix), gi(0)^V^,75g2(0)}|0) = -T,-^p;(g') iqz)% + ■■■ . (238) 



It is interesting to note that the above correlation functions receive no contributions from 
1~ states. The OPE results of T^p_^{q'^) have the following forms 



a 



-q^ {asG"^) W-n-tts 



7207r3 

as 
187rg2 

_,A,OPE/ 2^ 



q^\n— T^izrjT + + feg2)^) + ^7rr(gigi)(g2g2)) 



/i^ 367rg2 27g^ 



mi{qiqi) + m2{q2q2)^ - (rni{q2q2) + m2{qiqi)^ (^In-^ + ^ 



, (239) 



,p^ {q') = ■ {asG') - 



187rg2 



l{mi{qiqi) - m2{q2q2)^ +4^mi(g2g2) - m2{qiqi)^ ( In 



2 2 



. (240) 



_,o-,OPE^ 2 



367rg^ 



5(mi(gigi) -m2(g2g2)) +4(mi(g2g2) -"^2(gi'?i)) In 



-g2 7 



/i2 12 



+0 ■ {a,G') + ■ (gigi)(g2g2) + (0 + 0(«,')) (gigi)^ (g2g2)' • 



(241) 



50 



Consequently, we obtain the sum rules 
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Note that the calculation of Eq. (12361) is actually analogous to that of for the p p], where 
there is no 75 for the p. Neglecting the mass corrections in Eq. fl242p . our perturbative and 
gluon-condensate contributions agree with Eq. (C.15) in Ref. [9j, but the term of dimension-6 
is different from theirs, where the sign difference due to 75 has been considered. 



5.3 Results 

In the numerical analysis, we use /sp^, f^p_^, s^^^, Sf/'^, and the masses for axial- vector 
mesons, which have been obtained in the previous section, as inputs. We also adopt the 
parameters given in Appendix O The Borel window can thus be determined by means of 
that the contributions arising from the higher resonances, including the continuum, and from 
the term of the highest dimension in the OPE series are well under control. We find that 
the suitable Borel window is 1.5 GeV2 < M"^ < 2.5 GeV2 except that the Borel window is 
2.5 GeV2 < M2 < 3.5 GeV2 for the uYp sum rule. The reason that we have to choose a 



51 



higher Borel window for the ujYp_^ sum rule is because a lower Borel mass will lead to a slowly 
convergence at the quark-gluon level; for instance, at = 2.0 GeV^, the highest dimension 
(dimension=6) term still gives a large correction, ~ 33%, at the quark-gluon level. Instead, 
during 2.5 GeV^ < < 3.5 GeV^, the contribution arising from the highest dimension 
term is about 14% ~ 5%. On the other hand, because the average gluon momentum fraction 
in an axial- vector meson is not less than zero, therefore we should have ujYp_^ > —4, so that 
the region for < 1.6 GeV^ is strongly disfavored by the ujYp_^ sum rule. 

As we evaluate G-parity invariant parameters for l^Pi states, the corrections receiving 
from l^Pi states are negligible, and vice versa. Nevertheless, if calculating G-parity violating 
parameters for the Kia (Kib) state, the corrections originating from the Kib {Kia) state 
cannot be ignored. In the numerical study we use ^'^^"^ and ajj'^^^, given in Table HI 
as inputs. To obtain the relevant parameter in a sum rule, we replace the other twist-3 
parameters using the corresponding sum rules. For instance, in evaluating aYp^ , we substitute 
fYsp^ and /^ipj by the expressions given in Eqs. f ll63p and (12141) into Eq. (11651) . Note that for 
ujYp^, because the working Borel window is different from the others, we thus adopt directly 
the values of f^3p^ given in Table [5] in the study. 

We summarize the numerical results in Tables [5] and [6l where the theoretical errors are 
due to variation of all inputs and the predictive uncertainty within the Borel window. To 
illustrate the sum rule results, we plot the parameters for the ai(1260), i^iA as functions of 
the Borel mass squared in Figs. [Hand bi{1235),KiB in Fig. [5l where the central values of 
input parameters given in Tables El [3l HI and in Appendix [C] are used. For simplicity, we 
do not plot the results for fi,fs,hi and h^, which are analogous to the above ones. The 
f^3p^ sum rule results obtained from Eq. (I199p and from Eq. (I193p are depicted. Because 
the sum rule given in Eq. (11930 exhibits badly convergent behavior (see the discussion after 
Eq. (I195p ). we therefore choose to use the diagonal sum rule given in Eq. (I199p . where the 
lowest lying pseudoscalar contribution is substituted by using Eq. (I200p . Although the sum 
rule in Eq. (11990 cannot determine the sign of /^3p^, however there are two reasons that its 
sign should be negative. First, Eq. (I193p yields a negative f^3p^ at the large M^. Second, only 
a negative f^sp^ can result in a physical value of ui-p^ that should satisfy < ui-p^ < 16/3. 

Finally, it should be noted that the main contributions of some sum rules for parameters 
are due to the quark and quark-gluon condensates. These sum rules may suffer from the 
contributions of higher resonances at hadronic level and radiative corrections at quark-gluon 
level. The two effects may partially cancel each other and give about < 10% corrections to 
the numerical results for parameters. See also the discussions in subsection I4.4.3[ Here we 
do not include these effects in Tables [5] and El 
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Table 5: Parameters for twist-3 distribution amplitudes of l^Pi mesons at the scales /i = 1 
GeV and 2.2 GeV (shown in parentheses), where fy^p^, f^sp-^i and f^3p^ are in units of GeV^. 





ai(1260) 


/i 


/s 


KiA 


fV 

h,3Pi 


0.0055 ±0.0027 


0.0055 ±0.0027 


0.0054 ±0.0027 


0.0052 ±0.0027 


(0.0036 ±0.0018) 


(0.0036 ±0.0018) 


(0.0035 ±0.0018) 


(0.0034 ±0.0018) 




-2.9 ±0.9 


-2.8 ±0.9 


-3.0 ±1.1 


-3.1 ± 1.1 


(-2.9 ±0.9) 


(-2.8 ±0.9) 


(-3.0 ±1.0) 


(-3.1 ± 1.1) 










-0.13 ±0.16 
(-0.13 ±0.16) 




0.0022 ±0.0009 
(0.0012 ±0.0005) 


0.0022 ±0.0009 
(0.0012 ±0.0005) 


0.0028 ±0.0009 
(0.0015 ±0.0005) 


0.0026 ±0.0013 
(0.0014 ±0.0007) 


\ip^ 


— 


— 


— 


0.57 ±0.39 
(0.70 ±0.46) 


a^p^ 


— 


— 




2.4 ±2.0 
(2.4 ±2.0) 




-0.013 ±0.002 
(-0.009 ±0.002) 


-0.012 ±0.002 
(-0.009 ±0.002) 


-0.012 ±0.002 
(-0.009 ±0.002) 


-0.012 ±0.002 
(-0.009 ±0.002) 


Ujtp^ 


-3.7 ±0.4 


-3.4 ±0.4 


-3.2 ±0.6 


-3.4 ±0.6 


(-2.9 ±0.3) 


(-2.6 ±0.3) 


(-2.4 ±0.4) 


(-2.6 ±0.4) 


alp^ 








0.07 ±0.21 
(0.05 ±0.15) 



6 Models for LCDAs 

6.1 Two-parton LCDAs of twist-two 

We have calculated the first few Gegenbauer moments of leading-twist light-cone distribution 
amplitudes of l^Pi and l^Pi mesons in Sec. 14.41 using the QCD sum rule technique. The 
Gegenbauer moments of higher conformal spins may not be predictive in the QCD sum rule 
approach owing to the divergence of the OPE series for relevant correlation functions. There- 
fore the models for light-cone distribution amplitudes depend on the truncated conformal 
expansions with the reliable Gegenbauer moments. 

Here we take into account the approximate forms of twist-2 distributions for l^Pi mesons 
as follows: 



$11 (m) = Quu 
$^(m) = Quu 



l + 3a|e±a|-(5e'-l) 



ao^±3a^e±4 2(5^'-l) 



(245) 
(246) 
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Table 6: The same as Table [5] except for l^Pi mesons, where f^ip^^, f^ip-^, and /^ip^ are in 
units of GeV^. 





6i(1235) 


hi 


^8 




fV 


0.0052 ±0.0018 
(0.0030 ±0.0011) 


0.0046 ±0.0021 
(0.0027 ±0.0012) 


0.0045 ±0.0020 
(0.0027 ±0.0012) 


0.0049 ±0.0021 
(0.0029 ±0.0012) 




— 
— 


— 
— 


— 
— 


0.07±0.19 
(0.09 ±0.24) 




— 

— 


— 

— 


— 
— 


0.35 ±0.73 
(0.31 ±0.68) 


fA 


-0.0058 ±0.0023 


-0.0053 ±0.0023 


-0.0055 ±0.0023 


-0.0065 ±0.0029 


(-0.0036 ±0.0014) 


(-0.0033 ±0.0014) 


(-0.0035 ±0.0014) 


(-0.0041 ±0.0018) 




-1.5 ±0.4 


-1.9 ±0.6 


-3.5 ± 1.0 


-1.9 ±0.6 


(-1.4 ±0.3) 


(-1.7 ±0.4) 


(-2.9 ±0.8) 


(-1.7 ±0.4) 










-0.06 ±0.05 










\—\j.\jD ± U.U4j 


ft}Pi 


0.011 ±0.006 


0.012 ±0.006 


0.012 ±0.005 


0.012 ±0.005 


(0.006 ±0.003) 


(0.006 ±0.003) 


(0.006 ±0.003) 


(0.006 ±0.003) 










0.17±0.17 
(0.25 ±0.25) 










-0.71 ±0.53 
(-0.76 ±0.56) 



where ^ = 2u — l and a^, a^, are non-zero only for strange mesons. The above LCDAs of 
pure l^Pi states are normalized as the normalization conditions 



du(^\\{u) 



du^Au) 



(247) 



(248) 



In our convention, u is the momentum fraction carried by the qi quark in an axial-vector 
meson (and is therefore equivalent to the momentum fraction carried by the s quark in a 
strange meson). Note that in this paper the strange mesons that we discuss should contain 
an s quark, while for the LCDAs of strange mesons involving an s, the replacement u ^ 1 — u 
has to be made, namely ^ — ^. As for l^Pi mesons, we take the following approximation 
(see also the discussions given in Ref. 



$11 (m) = Quu 



$_l(m) = Quu 



.11 



+ ?>a\i + al-{he-l] 



1 ± 'ia^ { ± flj 



(249) 



(250) 
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Figure 4: Some relevant parameters in determinations of three-parton distri- 
bution amplitudes of twist-3 for the Kia (solid curve) and ai(1260) (dashed 
curve) as functions of the Borel mass squared, where the central values of input 
parameters have been used. The renormalization scale is set at yU ^ 1.4 GeV 
except that uYp^ is at /i ~ 1.7 GeV. For f^3p^, the lower two curves are derived 
from Eqs. fll99p and fl200p . and the upper two curves from Eq. fll93p . 
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Figure 5: Some relevant parameters in determinations of three-parton distri- 
bution amplitudes of twist-3 for the Kib (solid curve) and 6i(1235) (dashed 
curve) as functions of the Borel mass squared, where the central values of input 
parameters have been used. The renormalization scale is set at ~ 1.4 GeV. 
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where and are non-zero only for strange mesons, so that the LCDAs of pure l^Pi 



states are normahzed as the normahzation conditions 



du^\\{u) = al, (251) 



/o 



/ du^±{u) = 1 . (252) 
Jo 

Due to mixtures, we define the LCDAs of physical /ii(1170), /ii(1380), /i(1285), /i(1420) 
mesons in Appendix [D], and, on other hand, the LCDAs for i^'i(1270) and fri(1400) are 
consequently given by 

$^^(1270)^^) ^ fK,.rnK,, fK^n^K^n (253) 

fKi{1270)'m'Ki{1270) JXi (1270) "^i^Ti (1270) 

$f (^^°°)(n) = $;^^-^(^)cosgK- , ^^-^^^-^ ^f^^u)smeK, (254) 

Ji^i (1400) "^Xi (1400) Ji^i (1400) "^7^1(1400) 



^ji(i270)(^^ = fr^ 'l'f-^(M) sin^;^ + -^^i^$^is(m) cos^^, (255) 

•'ft:i(1270) JKi{1270) 



^Ki(i4oo)(^^ = ^(m) cos^/^ - -p^^^'l'^u) sinOK. (256) 

/a'i(1400) //i:i(1400) 

In Figs. E] and [71 we plot the twist-2 light-cone distribution amplitudes for l^Pi and 
l^Pi states, including results for the physical mesons, /ii(1170), /ii(1380), /i(1285), /i(1420), 
7^1(1270) and i^i(1400), at the scale /x = 1 GeV. ^_^^{u) and $||^^(m) are symmetric under 

M ^ 1 — M if neglecting SU(3) breaking effects, whereas $||^^(m) and $/'^(m) are antisym- 
metric. 

The contents of /ii(1170) and /ii(1380) are dominated by their uu and ss components, 
respectively, which are slightly different for 61 p-^ = 10° and 45°. However, the ss component 
of /ii(1170) and uu component of /ii(1380) become significant for 61 = 10°. Analogously, 
considering the real states /i(1285) and /i(1420), we find that only the ss fraction of /i(1285) 
is a little sensitive to the singlet-octet mixing angle 03 p^ as changing 03 p^ = 38° to be 50°. 
In particular, the ^± and $|| for /i(1285) (or for /i(1420)) are predominated by their uu 
(or ss) component. In Fig. Wi^), it is interesting to note that due to the mixture between 
KiA and Kip, where the axial- vector mesons contain an s quark and a light anti-quark q, 
we find that the q (or s) carries a larger momentum fraction for the i^i(1270) (or i^i(1400)) 
meson with respect to ^ = 45°. Nevertheless, if = —45°, the s (or q) instead carries a 
larger momentum fraction for the i^i(1270) (or i^i(1400)) meson. 

6.2 Two-parton LCDAs of twist-three 

Using the equations of motion allows one to rewrite the two-parton LCDAs of twist-three in 
terms of the leading-twist LCDAs and three-parton LCDAs of twist-3. Thus, substituting the 
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0.2 0.4 0.6 0.8 1. 0.2 0.4 0.6 0.8 1. 

u u 

Figure 6: Leading-twist light-cone distribution amplitudes, normalized at the 
scale /i =1 GeV, for l^Pi states, where the central values of Gegenbauer mo- 
ments given in Table H] are used, u {u = 1 — u) is the meson momentum fraction 
carried by the quark (antiquark). In (a) and (c), the solid, long-dashed, short- 
dashed and dot-dashed curves correspond to bi{1235), hi{singlet),hs{octet) 
and KiB, respectively. In (b) and (d), the solid [short-dashed] and long- 
dashed [dot-dashed] curves correspond to the mm [ss] contents of /;,i(1170) 
and /ii(1380), respectively, where Oip-^ = 10° (45°) have been used for heav- 
ier (lighter) curves. The definitions for the LCDAs of /ii(1170) and /ii(1380) 
have been given in Appendix [Dl 



twist-2 LCDAs specified by Eqs. f l245D . fl246D . fl249l) . and fl250l) . and three-parton LCDAs of 
twist-3 specified by Eqs. (H0l)- (H5l) into Eqs. ( l55l) . (1561) . ( |67j) . and ( |68l) . we get the approximate 
expressions in linear in quark masses (valid up to conformal spin 9/2): 



"S3,3Pi'^3Pi + y,3Pi y 



21 



e(5r - 3) 



9_ 



~ o '^1 ( o + + In M + In -u — -a^ 5_ (3^ + In -u — In -u 
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Figure 7: Leading-twist light-cone distribution amplitudes, normalized at the 
scale fi =1 GeV, for l^Pi states, and for ii'i(1270) and ii'i(1400) mesons, 
where the central values of Gegenbauer moments given in Table H] are used, u 
(m = 1 — m) is the meson momentum fraction carried by the quark (antiquark). 
In (a) and (c), the solid, long-dashed, short-dashed and dot-dashed curves 
correspond to ai(1260), fi{singlet) , fs{octet) and Kia, respectively. In (b) and 
(d), the solid [short-dashed] and long-dashed [dot-dashed] curves respectively 
correspond to the uu [ss] contents of /i(1285) and /i(1420), where 9 = 38° (50°) 
have been used for heavier (lighter) curves. The definitions for the LCDAs 
of /i(1285) and /i(1420) have been given in Appendix [Dl In (e) and (f), the 
solid (dashed) and dashed (sohd) curves correspond to i^ri(1270) and i^i(1400), 
respectively, for Ok = 45° {6k = -45°). 
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for pure ^Pi states. Note that to include the corrections consistently in linear in quark masses, 
in Eqs. f l257l) - fl260p and (12641) - fl267p the parameters with the "6ar" should be replaced by the 
corresponding ones in the massless quark limit. For the physical /;.i(1170), /;.i(1380), /i(1285) 
and /i(1420) mesons, their two-parton LCDAs of twist-3 are defined in Appendix [Dl The 
LCDAs for iri(1270) and iri(1400) are given by 

fx^A^K^A Ja,v)K^A a , fK^B^K 



^ia,v)Ki{1270) 



/Xi (1270) '^Xi (1270) 



gj_ sm + 



Ik^A^K: 



/Xi (1400) "^Xi (1400) 



lA {a,v)K^A n 

fi'_L COSt^x — 



IB gMK.B^^^Q^^^ (271) 

^7^'^)^-sin^^, (272) 



Ik 1 (1270) '^i^i (1270) 

Jk^b'^K^b {a,v)K 
/ftTi (1400) "^Xi (1400) 
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^(t,p)/..(1270) ^ fk^^^A j^MK^A ^.^Q^ ^ fks<.s J^(^,P)K,, ^^^^^^ ^273) 

J Xi (1270) "^ftTi (1270) f i^i (1270) "^Xi (1270) 

^(,p)K,(i4oo) ^ ^^^<i^^(,p)/.,. ^^^^^ _ /i.B^x,. f^(,,p)K,, ^274) 

^ (1400) ''^Ki (1400) ^ (1400) "^i^i (1400) 

Substituting the central values of parameters given in Tables Hj [5l and [6] into the above 
equations, in Figs. [HlfTT] we plot two-parton LCDAs of twist-3 at the scale fi = 1 GeV. 
The properties of twist-3 two-parton LCDAs are analogous to the cases of the leading-twist 
LCDAs that we have given previously. It is interesting to note again that the LCDAs, 
9±\9±\h^\ \ ^-iid h\^^\ for /i(1285) (or /i(1420)) are dominated by the uu (or ss) content 
which is insensitive to singlet-octet mixing angle dap-^ in the range 38° < 9ip^ < 50°. More- 
over, gf^'', and /ly^^ for /ii(1170) (or /;,i(1380)) are dominated by the uu (or ss) 
content for 9ip^ = 45° but the ss (or uu) content becomes significant for 9ip^ = 10°. 

6.3 Three-parton LCDAs of twist-three 

The approximate three-parton LCDAs of twist-3 are given in Eqs. ( l40l) -f H5l) and the relevant 
parameters are summarized in Tables O and O For completeness and simplicity, we plot the 
LCDAs for the ai(1260) and 6i(1235) mesons in Fig. [12] to illustrate their behaviors. 

7 Summary 

The light-cone distribution amplitudes specified by the collinear twist can be expanded in 
terms of the series of the so-called conformal partial waves. Each partial wave is characterized 
by a specific conformal spin. For each conformal spin, the dependence of the distribution 
amplitudes on the transverse coordinates is governed by the renormalization group equation, 
while the dependence on the longitudinal coordinates is involved in "spherical harmonics" 
of the SL{2,R) group. 

We have presented a detailed study of twist-2 and twist-3 light-cone distribution am- 
plitudes of axial-vector mesons, based on the QCD conformal partial wave expansion. The 
equations of motion allow us to obtain the relations among the twist-two and twist-three 
light-cone distribution amplitudes so that we can use a minimal number of independent 
nonperturbative parameters to describe the distribution amplitudes. The conformal partial 
wave related by equations of motion should correspond to the same conformal spin since 
equations of motion in the QCD perturbative theory respect all symmetries given in the 
classical level. 

Our main results are as follows: 
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Figure 8: Twist-3 light-cone distribution amplitudes g^f^ at the scale /x =1 GeV, 
where the central values of parameters given in Tables HI O and E] are used. In 
(e), where 6k = 45°, and in (f), where 6k = —45°, the solid and dashed curves 
correspond to -ft'i(1270) and ii'i(1400), respectively. Others are the same as 
Figs. Eland El 



• In subsections 14.21 and 14.31 we have shown the sum rule results for the axial- vector 
(tensor) decay constants of l^Pi (1^-Pi) axial- vector mesons, where we have updated 
the values for l^Pi states obtained in Ref. |16j . 

• Using Gell-Mann-Okubo mass formula, we have obtained the mixing angles for the 
/s (octet) and /i (singlet) of l^Pi states to be 6*3 ~ 38°, and for (octet) and hi 
(singlet) of l^Pi states to be 6ip^ ~ 10°. Thus the decay constants and light-cone 
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distribution amplitudes for these states are determined. 

• The sum rules for the first few Gegenbauer moments of the leading-twist light-cone 
distribution amplitudes together with their numerical results have been given in sub- 
section H31 where the SU(3) breaking effects relevant to the Kia and Kib states are 
included. The results for l^Pi states, hi,h^ and for G-parity violating Gegenbauer 
moments are new, while the results of G-parity invariant Gegenbauer moments for hi 
and KiB (which are l^Pi states) are updated. In the sum rules, because the G-parity 
violating Gegenbauer moments of l^Pi and l^Pi are always mixed togther, we thus 
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Figure 10: Twist-3 light-cone distribution amplitudes /ly at the scale 
fi =1 GeV. Others are the same as Fig. [HI 



add a reasonable constraint 

-,KiA 



° =1.0 ±0.3 



a. 



IS 



which is a good approximation for G-invariant Gegenbauer moments, to obtain quali- 
tative estimates. See the detailed discussions in subsection I4.4.3[ 

In Sec. O using the QCD sum rules, the relevant G-parity invariant and violating 
parameters for expanding the three-parton distribution amplitudes of twist-3 in terms 
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Figure 11: Twist-3 light-cone distribution amplitudes Iv^ at the scale 
\i=\ GeV. Others are the same as Fig. [81 



of conformal partial waves with conformal spin up to 9/2 have been evaluated, where 
the SU(3) corrections have been contained. To determined not only the magnitudes 
but also the relative signs for the parameters, one of the interpolating currents in the 
two-point correlation functions is chosen to be the local axial-vector (or pseudo-tensor) 
current in calculating the parameters for the l^Pi (or l^Pi) state. All the results are 
new. We have checked all the calculations very carefully since nobody did these before. 
Only the sum rule calculation for j^\p^ is very similar to that for in SU(3) limit 
(See the discussions after Eq. (I244p ). On the other hand, the resulting f^ip^ sum 
rule is not reliable since the calculated OPE series is not well convergent. We further 
resort to the diagonal sum rule for ftp\ and this calculation can be found in Ref. |20j . 
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Figure 12: Twist-3 three-parton light-cone distribution amplitudes for the 
ai(1260) and 6i(1235) mesons at the scale fi =1 GeV: (a) / = ^for ai(1260) or 
/ = V,r for 6i(1235), (b) / = V for ai(1260), (c) / = T for ai(1260), and (d) 
/ = ^ for 61(1235), where we have taken the central values of uj^_^,uj^_^, and 0;^, 
given in Tables [5] and El as inputs. Here and are momentum fractions 
carried by the quark and anti-quark in an axial-vector meson, respectively. 
The gluon momentum fraction is substituted as = 1 — — ag^. 

where the sum rule is used to study the coupling of the kaon. The diagonal sum rule 
is stable and the sign for f^-p^ can be determined in an indirect way (See Sec. l5.1.3l for 
more discussions). It should be noted that as G-parity violating parameters for l^Pi 
(l^Pi) states are computed, the corrections receiving from l^Pi (1^-Pi) states have to 
be considered. 

• Adapting the EOM formulas derived in Ref. [9] for vector mesons to the present case, 
the two-parton LCDAs of twist-three can be written in terms of the leading-twist 
LCDAs and three-parton LCDAs of twist-3. The detailed results for twist-three three- 
parton LCDAs are shown in Sec. 16.21 In the SU(3) limit, the detailed symmetric 
properties of LCDAs of axial-vector mesons, as compared with that of the vector 
mesons, are summarized in Table [H 
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• Using the conformal partial expansion, we presented the models for light-cone distri- 
bution amplitudes, containing contributions with conformal spin up to 9/2, in Sec. [61 

• We have considered the strange quark mass corrections to distribution amplitudes for 
the strange axial-vector mesons, /('i(1270) and i^'i(1400). It is interesting to note that 
^^f^{u), $11^^ (u), (7^^ g^^^ , h\^l'^ have significant antisymmetric behaviors, 
which should be phenomenologically attractive. 

It should be noted that because corrections due the higher resonance and radiative correction 
in OPE may partially cancel each other for non-diagonal sum rules, it is estimated that the 
resultant errors for parameters of LCDAs may be ~ 10%. On the other hand, because we 
do not calculate the radiative corrections to the perturbative term in diagonal sum rules, 
these corrections may also lead to ~ 10% errors for results. However, since we do not do 
the qualitative calculations about these effects, we thus do not include this possible error in 
Tables H El and El 

Recently, Belle has measured i^f (1270)7 s-^icl given an upper bound on 

ii'f (1400)7 [H]- Interestingly, the recent calculations [lH ^ of adopting LCSR (light-cone 
sum rule) form factors [15] gave too small predictions for B{B^ — > i^f (1270)7) compared 
with the data. Since the physical states -ft'i(1270) and -ft'i(1400) are the mixture of Kia 
and KiB which are respectively the pure l^Pi and l^Pi states, the light-cone distribution 
amplitudes of Kia and Kib are relevant to the results oi B ^ 7^1(1270) and i^i(1400) 
transition form factors. It is known that for Kib-, ^\\ is antisymmetric, while $_l is symmetric 
in the SU(3) limit due to G-parity. Nevertheless, for Kia-, $|| becomes symmetric, while $_l 
is antisymmetric. The above properties were not correctly studied in the literature. Some 
related researches will be published elsewhere [T7] . 
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Appendices 



A Spin Projections and Collinear Twist 

Considering infinitesimal rotation — > = + e^^x" in the four dimension, tlie general 
field $(x) transforms as <l>'(x') = [1 - e^^(x^9'' - x^df" - S'^'')]$(x) = [1 + ie f,^M (x), 
where S^'^ is called the generator of spin rotations of the field $. For scalar, quark and gluon 
fields, we have 

S^Xx) = 0, = ^a^""^, WA"" = g'^^'A^ - g^'^'A" , (A.l) 

respectively, where A^ = T'^A'^^. We can apply the following spin projections 

n = ^7-7+, P- = ^7+7-, P+ + P- = l (A.2) 

on a quark field ifj to project its typical spin component, so that a quark field with a given 
spin component s on the moving direction can be measured to be 

S+_P±^ = s±P±^ = ±^P±^ . (A.3) 



Recalling that the canonical dimensioro of a quark field is ^ = 3/2, we can therefore de- 
compose an arbitrary light-ray quark-antiquark current into different {collinear) twists (= 
dimension — spin projection on the moving direction) components: 

V^2(x2-)r^l(Xi_) = ^2(x2-)P-rP+^i(Xi_) 

^ V ' 

twist-2 



V ■ 

twist-3 

+ ^2(x2-)P+rP_^i(xi_), (A.4) 
twist-4 

where F stands for a generic Dirac matrix structure. Note that, taking '?/'2(x2_)P-FP_^.^/'i(xi_) 
as an example, its {collinear) twist exactly equals to t = ii + £2 — si — S2 = 2, while its 
conformal spin is j = {£1 + £2 + Si + S2 + n) /2 = ji + j2 + n with n = 0, 1, 2, ■ ■ ■; in other 
words, j is not a fixed value for a LCDA defined by a non-local composite field as seen in 
Eq. ( |T3l) . On the other hand, for a gluon field with the canonical dimension £ = 2, one can 
find that 







= 1 






-G_i_± 


= 


• G_i_± 








1 ■ G_ 



S_| — ■ G-i , 

(A.5) 



10 



Here we do not distinguish the canonical dimension ^'^^^ and scahng dimension 



69 



where, for instance, G+j. = Gfj_i,n^g'^ and g'^ = g^" — n'^n°' — n^n^ . In analogy to the 
previous discussion, we can decompose an arbitrary hght-ray quark-gluon-antiquark current 
into currents with different twists: 

i)2{x2-)'^9sG^y{x^-)i)i{xi^) = i)2{x2-)P-TgsGap{x'i^)P+i)i{xi^){g°^n'^n'' + n'^n^'g'^^) 

twist-3 

+ ^2{x2-)P+Tg,G^p{x^_)P+'iP^{x^.){gTnPn'' + n'^n^gi') 

twist-4 

+ i^2MP-TgsG^p{x:,^)P^iPi{xr^){gTn^n' + n^n'^^^f ) 

^ v ' 

twist-4 

+ ^2(x2-)P-r(7,G,^(x3-)P+^i(xi^)((77(7f + n^n'^nV + n^n'^nV) 

^ V ' 

twist-4 

+ C(twist-5,6,7) . (A.6) 

B Operator Identities 

The operator identities, which are used to obtain the integral equations, Eqs. (IS^ . 
dMD, and dSSD, follows: 

d _ 

q2ix)^^qi{-x) = / dt- — q2{tx) /kqi{-tx) 
Jo ox^, 

dtt / dv q2{tx)gsG ^y{vx)x'' /c'-^^qi{—tx) 
J~t 

i dt dvv q2{tx)gsG^u{vx)x'^ /tqi{—tx) 
Jo J-t 



-if^i^vcn / dttx^'d" [q2{tx)'^'^'^5qi{-tx 
Jo 



+ {rriq^ - mq^)x'' dtt q2{tx)(Jy^qi{-tx) , (B.l) 







d _ 

q2{x)l^,imi{-x) = J dt^q2{tx) /k-f5qi{-tx) 



1 /■* 

dtt / dv q2{tx)gsG fj_y{vx)x'' ^qi{—tx) 
J-t 

i dt dvv q2{tx)gsGfj_u{vx)x'^ /t'y5qi{—tx) 
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Jo 



+ {ruq^ + mq^)x'' dttq2itx)a^i,-f5qi{-tx), (B.2) 
Jo 



d _ 

-Q^q2{x)(Tt,u'^m{-x) = -iduq2{x)'j5qi{-x) + J dvq2{x)x°'gsGau{vx)-f5qi{-x) 



-i J dvvq2{x)XpgsGP^{vx)ai,^-f5qi{-x) 

- {nig^ + mq^)q2{x)-f^-f5qi{-x) , (B.3) 

^2(a;)75gi(-a;) - g2(0)75gi(0) = dt dvq2{tx)x°'aai3X^gsG^'^{vx)-f^qi{-tx) 

Jo J-t 

+i / dtd"^ [q2{tx)(Tci3X^'^^qi{-tx)] 
Jo 

+i{mq2 - mq^) / dtq2{tx) ^-f^qi{-tx), (B.4) 

^0 



where we have adopted the following notation to stand for the total derivative: 



d 

{q2{x)Tqi{-x)] = — {q2{x + y)[x + y, -x + y]Tqi{-x + y)} 

oy^ 



(B.5) 



y^O 



Eqs. (IB.ip . (IB.2p . and (IB.3P have been obtained in Refs. [HIE], whereas Eq. (IB.4p is new 



C Input parameters 

The theoretical input parameters, used in our analysis, together with their respective ranges 
of uncertainty are summarized here. We take into account ^^(l GeV) = 0.497 ± 0.005, 
corresponding to the world average as{mz) = 0.1176±0.0020 [28], and the following relevant 
parameters at the scale /i = 1 GeV [i6| Il0l[20t [39] : 



{asC^^G^f"') = (0.474 ±0.120) GeVV(4vr) 
{uu) = {dd) = -(0.24 ± 0.010)3 GeV^ , 
{mu + md)/2 = (5±2) MeV , 
{QsUaGu) = (gJaGd) = -(0.8 ± 0.1){uu), (gsSaGs) = (0.8 ± 0.1)(^,u(tGu), 
f^pg = (0.0045 ± 0.0015) GeV^ , 
= 0.8 GeV^ , 



{ss) = (0.8±0.1)(mm) , 
m, = (140 ±20) MeV^ 



-,/8 

'0 



(C.6) 



Sn ^ St," ~ s^^ = 1.1 GeV^ 



,K* 



-,P _ 



= 0.04 ±0.03, 
0.2 ±0.1, 
= 0.13 ±0.08. 



= (185 ± 10) MeV 
aa'"^ = 0.0 ±0.1, 
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where the scale-dependence of operators is given by 

{gsQCT ■ Gq){Q) = {gsQa ■ Gq){^) 



{asG'){Q) = {asG^){lj,), (C.7) 

with b = {llNc — 2nf)/3. As described in Eq. fl86l) . we adopt the vacuum saturation approx- 
imation for describing the four-quark condensates, i.e., 

(OlgT.A'^ggr.A'^glO) = --L_Tr(r,rOTr(A'^A")(gg)2 . (C.8) 

Performing the analysis in analogy to that given in Ref. and using the results in Refs. [ITf 
ill SI] and in Eqs. ([WD, ([HH]), ([USD, (12®, ([209D, ([HQD, ([221, dZSSD, ([226D, we obtain 
the LO scale-dependence of the twist-3 parameters with the light quark mass corrections as 
below 

[G-parity invariant components]: 

/3^3P,(Q) = L^^^'fl^pM)^ r+ = -Ic^ + 3Gg, (C.9) 



V _|_ 28 
Pi 3 •'3 
y , ,V 28 fA 

T^3,3Pi 



3,3Pi'^3pj 



f: 



3,3Pi^3 



V _|_ 28 f A 
Pi 3 •'S.^Pi 

V 28 f A 
^3 



Pi 



T^3,3Pi 



^Gp + AGg |Cp — |C| 



2^ 2^ 8^ , 7^ 



(C.IO) 



^ = + Gg, 



(c.ii: 



/3t3p,c^3^P,(Q) = ^'^"''/V3t3p,c^3^p,(/i) - ^(l^/^ -L^'^'/^jK, +m,J(^)/3p, 



85 

r^3- = + -Cg, 

6 3 ' 



(C.12) 
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[G-parity violating components]: 

fs^pXp^iQ) = L^^^'ft^P^pM^ ^2 = -Icp + 3Cg, (C.13) 

28 fV _ fA \ ^ / 28 fV ^r^ _ fA ^A \ 

t/3,3Pi'^3Pi + fsSPi'^sp^ j ^ \ ■3"/3,3Pi<^3p^ + fsap^CTsp^ J ^ ' 

r+ ^ ( 3*^^ ^ ~ ^'^^ 1, (C.14) 

r^3+ ^^Cp + Cg, (C.15) 



for l^Pi states, and 
[G-parity invariant components]: 

2 

29 



~795 -LW{36) + i,-56o/(96) j Ki +^?2)(Ai)/iPi(A^), (C.17) 

/3^P,(Q)=i^''^"/V3tip,(/^), (C.18) 

[G-parity violating components] : 

fhp^'(pM = L^"'"fhpAM + - L^^/'')K, -m,J(/.)/.-^p^(/.),(C.19) 

{■V ,2&fA ^A \ ( fV „V _,28fA ^A 

^3,iPi'^iPi ^ 3 ^a.iPi'^iPi 1 _ T-rj/fe ( is.iPi'^iPi ^ 3 ^3,iPi iPi 

f V _28fA A I - I _28fA A 

/3,lPi^lPl 3 /3,lPi^lPi / Q \/3,lPi^lPi 3 h^Pi'^^Pi ; ^ 

7080 V L^^/^^'') - L^^^/(^^) j ("^91 "^q2)(/^)/ipi(/^)) 

(C.20) 

/3tiP,Afp,(g) = i^^"^^/V3tiP,Afp,(//), (C.21) 

/3tip,<^iPi(Q) = "it-pApM)^ (C-22) 
for l^Pi states, where Cp = (A^^^ _ i)/(2Ar^), Cg = A^c and L = as{Q)/as{fj). 
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D The definitions of LCD As of physical fi and hi mesons 

For /i(1285),/i(1420),/ii(1170) and /ii(1380) mesons, denoted by B in the following, we 
define the chiral-even LCDAs to be (with Q = u,d, or s) 

/■I f *{^) 



and the chiral-odd LCDAs to be 



+ 



{pzf 



(5(P,A)|g(i/)75g(x)|0) = /^'^m|(e*Wz) / ^/m e^("^'^+^P-) , (D.26) 



where the distribution amplitudes are subject to the choices of the normalization constants 
with ii) D = u for /i(1285), /ii(1170), and (ii) D = s for /i (1420), /ii (1380). Here 
we use the conventions for the decay constants that /^("285) ~ //i{i285)' //^'{i420) ~ //i(i285)' 
/hi(ii70) = /mii70)(1 GeV), and A'^(i38o) = ftlm^^ relevant decay constants 

can be found in Eqs. (I9l])-(l9l]) and ffTO^ -flM. 
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